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(54) Method for determining a partition in order to insert a watermark, and associated insertion 
and decoding methods 

(57) A method of determining a partition into distinct 
regions of a set of coefficients representative of at least 
a part of an image in order to insert a watermark com- 
prises the following steps: 

generation (E501 ) of a centred pseudo-random se- 
quence (w) equal in size to the cardinal number (N) 
of the said set of coefficients, formed from centred 
pseudo-random sub-sequences; 

- modulation (E502) of the said set of coefficients by 
the said centred pseudo-random sequence (w) in 
order to insert the same information bit on the said 
set of coefficients; and 

- checking (E504-E508), over each region, of a crite- 
rion of detectability of the information bit inserted by 
modulation. 

Use especially for determining an adaptalive parti- 
tion as a function of a criterion of detectability of an in- 
formation bit inserted on each region. 
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Description 



roooil The present invention concerns a method of determining a partitioning ol a digital image for inserting i a wa- 
SSSLS^m-hod of inserting a watermark, a method of processing an image and an assoaated method of 

5 ^rSo^Z^Tc^s a device for determining a partitioning of a digital ,mage for inserting a watermarking 
ESS a ^oMnserting a watermark, a devbe for processing an image and a dev.ce for decod,ng a watermark.ng 

[W03] The present invention lies in general terms in the technical field of watermarking of digital images, and more 

10 g£ffl£^ digital data makes it possib.e to protect these data, for example by associating copyright info, 

m "in"" genera, principe. watermarking consists of inserting an indelible watermark in digital data, which can be 

ss x v^^^ a ^ e ^ p— - - — — ' has been 
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inserted. 



[00071 This inserted watermark must consequently be both imperceptible, robust to certain distort.ons applied to the 

^TZ^^s^^ of inserting a watermarking signal in a digital image consists of using a 
ffnSS I which at least one subset of coefficients representing the digital image .s modulated ac- 

^r*d^ 

ST! { V < i < PI - watermark of size P <N,a pseudo-random signal with known distribution and average of zero, 
the linear insertion formula is : 



XV = Xj + b otj Wj, with 1 < j < P, 
in which {Xj, 1 < j < P) is a subset of the set of coefficients X. b isan information bit, and a, isa weighting coefficient. 

30 Th^iec^Xtatermark then consists in detecting whether or not the pseudo-random sequence w 

Z Sen insertea into a set of coefficients. This detection is carried out without using the or.g.nal .mage and can be 
based on L sStda d sed statistical test which makes it possib.e to calculate a probability of correct ^ectoa 
15)11] Such an insertion technique makes it possib.e, by insertion of a watermark, to .nsert a s.ng.e .nformat.on M 

s i S for example, the name or the address of the owner or of the author of the .mage 

loSSe the previously described insertion method as many t,mes as there are .nformat.on b,ts to be .nserted. 

nosed on one another, which could disturb the detection or cause troublesome v.sual effects. 

[MIS] ItTs a matter, consequently, of choosing a partition of the coefficients representee of the d.g.ta. .mage .nto 

45 ^i:^^"^*'^* inserting a watermarking signal o, given size by spectrum 

' The b aim of the present invention is to propose a method of determining a partitioning of a digital image which 
s^dapted to mTimage, ensuring detectabi.ity of the bits inserted for a technique of insert.on by modu.at.oa 
ss Wit] The present invention thus relates to a method of determ.ning a partitioning .nto a*^"*™^*^ 
mage in order to insert a watermarking signal by modulation of coefficients of said reg^ns represent^ a d grtat mag*. 
SZq in accordance with the invention, this partitioning into regions is effected by an adapts part.ton.ng as a 
function of a criterion of detectability of an information bit inserted m each reg.on. 
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sAss - ~i» « — - — -* 

eients. lire more reliable is the deleetrov calculation k usee lor ejecting me deleclioa the response 

[0023] II oar. In lacl be shown thai, when a 00 relation cs Kraal .0 mtemarkin g siajral 

!rooe=;:=;:'~ 
Seb=;^=r^== 

10028] This .s because. Ihe g.earer the moowatsav ^ u *' *2,^° c „». nmHW 'he greater 

"Jore reliable I. are d.reolion In that . ? fe Set m.^aeanarked ln»ge. 

signal and that it decreases with th.s length. ensuring imperceptibility of the watermarking signal 

K„«m^^^ 

5S panalaning mue pedoryr*d „eo ensure, MM* o. Ihe nrodu.,.*. e.eored on each o, ,h. r.grons 

obtained by the partitioning. . harar , P ri S tic of the invention, the detectability criterion is calculated using, 

r0O331 According to an advantageous characteristic of the invemi , = whefe k . |s a 

„ v is «iual 10 said nraxinrrinr value ol Ihe r^ghUrrg, coefcien , Clarion by Ihe inrage «gnel MH. 

45 a prior step of comparing the s.ze of each region of the P ar ™'"9 ed jn sajd regj0 n. 

,0039] According .0 one panlcul-^ P«euc^b«»n»^M «W««bWy s^pUud. I. cafcuWJd 

L»oi«,. S end »,e wel.rnrerking eigne, o, .he '""^ ° ^*£,^ ™s. in Prop. ICIP. pages 215-2.8. Sep. 
mo article enttled "A method lor signarore casing a, d.g.ral mages by 
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delegability of the inserted watermarking signal in less than or equal to the cardinal number of sa.d region. 
[0042] According to a preferred embodiment of the invention, the method comprises the following steps 

- generation of a centred pseudo-random sequence equal in size to the cardinal number of a set of coefficients 
s representative of at least a part of the image, formed from centred pseudo-random sub-sequences: 

- modulation of the said set of coefficients by the said centred pseudo-random sequence in order to insert the same 
information bit on the said set of coefficients; and 

- checking, over each region, of a criterion of detectability of the information bit inserted by modulation. 

io . [0043] It is thus possible to achieve a partition of a set of coefficients into distinct regions, by carrying out an adaptive 
oartrtioningonthebasisof a criterion of detectability of an information bit inserted on each region. 
[0044] Such an adaptive partitioning exhibits the advantage of adapting to the content of the image signal and of 
guaranteeing the detectability of the watermark inserted on the different regions determined by this partition, in contrast 

to the choice of a random partitioning.' ...... rf ^ K .. 

is [0045] In practice, for each region of the partition, a criterion of detectability of an information b.t inserted by modu- 
lation of the coefficients of the said region is checked. This criterion may consist in the comparison of a result of a 
standardised statistical test with a threshold value corresponding to a probability of detection fixed in advance. 
[0046] By virtue of the method in accordance with the invention, it is possible to generate a centred pseudo-random 
sequence on the basis ot the merging of centred pseudo-random sequences of smaller sizes, and to modulate all the 

20 available coefficients in order to achieve the insertion of a watermark. 

[0047] By modulating all the coefficients with the same information bit and a centred pseudo-random sequence which 
is broken down into pseudo-random sub-sequences which are also centred on each region of the partition produced 
over the set of coefficients, the steps of generating a pseudo-random sequence and of modulation do not have to be 
reiterated at each new region tested. 

2S [0048] According to one preferred characteristic of the invention, the method of determining a partition comprises a 
single step of applying a distortion to the said set of modulated coefficients, before the step of checking of a detectability 

[0049] " It is known to apply a predetermined distortion to a set of modulated coefficients and to carry out the check 
on the detectability criterion after this pre<Jistortion, in such a way that it is possible to guarantee the detection of the 
30 watermark inserted on a partition of the coefficients for a predetermined distortion 

[0050] By virtue of the invention, a pre-distortion can be applied to the set of modulated coefficients, in a unique way. 
and does not need to be reiterated at each new region tested in the partition 

rO051] In one particularly practical embodiment of the invention, the set of coefficients is a set of spectral coefficients 
originating from a spectral transformation of a digital image, and the step of applying a distortion includes the following 
35 sub-steps. 

- inverse spectral transformation of the said set of modulated coefficients in order to obtain a watermarked image; 

- application of a distortion to the said watermarked image; and 

- spectral transformation of the watermarked image in order to obtain a set of modulated coefficients after distortion. 

40 

[0052] The saving obtained in terms of complexity of calculation is all greater in this case whenever the application 
of a distortion requires several steps ot processing of the set of coefficients, and particularly the application of an 
inverse spectral transformation in order to obtain the watermarked image and to apply a distortion, such as compres- 
sion/decompression to it, and then the application of a forward spectral transformation of the image again so as to 
4S check the criterion of detectability in the spectral domain. 

[0053] According to another particularly practical characteristic of the invention, the centred pseudo-random se- 
quence is formed from pseudo-random sub-sequences of identical size. 

[0054] Such a pseudo-random sequence is particularly well adapted to the applicalion of a bottom-up or top-down 
iterative partitioning of the quadtree type which makes it possible iteratively to combine regions of identical sizes or 
so iteratively to divide regions into sub-regions ot identical sizes. 

[0055] For preference, the pseudo-random sub-sequences have a size greater than or equal to a minimum number 
of coefficients forming a region of the said partition of a size which is statistically significant for validating the sa.d 
criterion of detectability over the said region. 

[0056] It is thus possible theoretically to form a maximum partition of the set of coefficients with the same unique 
ss centred pseudo-random sequence, each pseudo-random sub-sequence being centred on the regions with the smallest 
size which is statistically admissible for inserting a watermark by modulation into the set of coefficients. 
[0057] According to one practical embodiment, which is inexpensive in terms of calculations, at the generation step_ 
each centred pseudo-random sub-sequence is created by generating half of the pseudo-random samples of the said 
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sub-sequence via a pseudo-random number generator and by symmetrising the said pseudo-random samples gener- 

[0056] For preference, in this embodiment, at the generation step, each centred pseudo-random sub-sequence is 
created by moreover permutating the said pseudo-random samples using a secret key. 

[0059] This permutation operation makes it possible better to preserve the intrinsic random-sequence nature of a 
sub-sequence created by symmetrising pseudo-random samples, and not to diminish the imperceptibly of the inserted 

rooeorAccording to another practical embodiment of the invention, at the generation step, each centred pseudo- 
random sub-sequence is created by generating pseudo-random random samples via a pseudo-random number gen- 
erator and by redistributing over each sample the sum of the said pseudo-random samples generated. 
[00611 According to a preferred characteristic of Ihe invention, a partitioning by iterative d.v,sion of the digital .mage 
!s performed and. for each region al a rank of said partitioning, a partitioning of a directly higher rank of said region is 
validated if and only A there exist at least two subregions of said region for which said detectability criterion is validated^ 
[0062] Thus, when it is attempted to divide a region into subregions, this partitioning of a higher rank is validated 
only if the number of regions in which a watermarking signal can actually be inserted whilst complying with the detecta- 
bility criterion is actually increased rt .,„ ta 
[0063] Otherwise, it is preferable to keep the region before division, of greater size, in order to increase the detecta- 
bility of the watermarking signal. ar1it . inn , h -wa| llP 
[00641 According to a preferred embodiment of the invention, for each region at a rank of said partitioning^ value 
of a capacity ol the digital image is incremenled when the partitioning of directly higher rank is not validated, the said 
capacity being equal to the cardinal number of the distinct regions obtained by the said partitioning on which the de- 
tectabilily criterion is validated. Ih .»,-h.rni 
[0065] It is Ihus possible to calculate the capacity of the digital image which thus corresponds to the number of 
information bits which can be inserted in the digital image. 

r00661 By virtue of the adaptive partitioning performed by the invention, by reiterating the partitioning of each region 
as long as the detectability criterion is complied with, there is obtained the maximum capacity of the digital image which 
corresponds to the watermarking capacity of this image by the technique of insertion by modulation. This watermarking 
capacity can be defined as the maximum number of information bits which can be inserted imperceptibly in the .mage 
with a guarantee of detectability on decoding. 

[0067] According to another embodiment of the invention, for each region at a rank of said partitioning, the value of 
a capacity is incremented when the partitioning of directly higher rank is not validated, said value of the capacity is 
compared with a fixed capacity value and the partitioning is reiterated for a region having a detectability amplitude 
which is the highest amongst all the other regions to be processed when said capacity value is less than said fixed 

[0068] ty ^partitioning method makes it possible, in this embodiment, to find the best possible partitioning in terms 
of detectability when a predetermined number of information bits have to be inserted, equal to said fixed capac.ly value^ 
[00691 By selecting, during the partitioning, the regions for which the detection amplitude is the highest, the present 
invention makes it possible to determine the partitioning which, a priori, guarantees the best detectability on decoding 
of the watermarking signals making it possible to insert a predetermined number of information bits. 
[00701 According to another embodiment of the invention, a partitioning is effected by iterative merging of the digital 
image and, for each region of said partitioning for which said detectability criterion is not validated, said region is merged 
with another region of said partitioning 

[0071] Preferably, in this embodiment, said other region of the partitioning is if possible a region for which said de- 
tectability criterion is not validated. H«««*flhiiHu 
[0072] It is thus possible to increase the capacity of the digital image by joining the regions in which the detectability 

criterion is not satisfied. 

[0073] According to another aspect of the invention, a method of inserting, into a digital image, a watermark com- 
prising one or more information bits, comprises a method of determining a partition in accordance with the first aspect 
of the invention, and a step of inserting the information bit or bits by modulation of the coefficients of one or more 
so reqions of the said partition. . 
F0074] This insertion method exhibits advantages similar to those of the method of determining a partition which it 
implements, and in particular makes it possible to insert a watermark in an adaptive and reliable way as regards its 

[SnS^Jhe present invention also concerns a method of decoding a watermarking signal in a digital image, inserted 
ss by modulation of coefficients representing said image in distinct regions forming a partitioning of sa.d digital image 

r0076] In accordance with the invention, this decoding method includes a step of determining a part.t.on.ng into 
distinct regions of the digital image to be decoded by effecting an adaptive partitioning as a function of a criterion of 
detectability of an information bit obtained by demodulation of the coefficients on each region. 
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[0077] Thus this decoding method makes it possible to find, at the decoder, the adaptive partitioning perlormed for 
inserting a watermarking signal. 

[0078] This is because, when the partitioning has been perlormed in an adaptive fashion as described above by a 
method of determining a partitioning in accordance with the first aspect of the invention, this partitioning, which depends 

5 on the content of the image, is not known to the decoder. 

[0079] The decoding method according to the invention then consists of testing all the acceptable partitionings and 
verifying the detection of a watermarking signal on each of the regions associated with this partitioning 
[0080] According to a preferred characteristic of the invention, the decoding method includes a prior step of comparing 
the size of each region of the partitioning with a minimum size corresponding to the minimum size of a statistically 

10 significant sample for the detection of a watermarking signal inserted in said region. 

[0081] According to a practical embodiment of the decoding method according to the invention, for each region of 
the partitioning, a detection amplitude is calculated from said detectability criterion and said detectability criterion is 
validated by comparing said detection amplitude with a threshold value predetermined for decoding. 
[0082] The decoding method thus makes it possible to find the adaptive partitioning effected when the watermarking 

15 signal is inserted. 

[0083] According to advantageous characteristic of this embodiment in the invention, the threshold value predeter- 
mined for decoding is less than said predetermined threshold value used during the method of determining a partition- 
ing. 

[0084] This is because, at decoding, the detection test is performed on coefficients which have a priori been modu- 
lo laled and possibly made noisy. These coefficients are different from the original coefficients, and the result of the 

detectability test at decoding will be different from that perlormed in order to determine the partitioning. By choosing 

different thresholds on decoding and coding, it is possible to take into account the tendency for variation in detection 

amplitude due to an increased variance in the coefficients which have been modulated and possibly made noisy. 

[0085] Correlatively, the present invention also concerns a device for determining a partitioning into distinct regions 
25 of a digital image for inserting a watermarking signal by modulating coefficients of said regions, comprising partitioning 

means adapted to effect an adaptative partitioning as a function of a criterion of detectability of an information bit 

inserted on each region. 

[0086] According to a preferred embodiment, this device for determining a partition includes: 

30 - means for generating a centred pseudo-random sequence equal in size to the cardinal number of the said set of 
coefficients, formed by centred pseudo-random sub-sequences; 

means for modulation of a set of coefficients by the said centred pseudo-random sequence in order to insert the 
same information bit on the said set of coefficients; and 

means for checking, over each region, a criterion of detectability of the information bit inserted by modulation. 

35 

[0087] The present invention also envisages a device for inserting into a digital image a watermark comprising one 
or more information bits, comprising a device for determining a partition in accordance with the invention, and means 
for inserting the information bit or bits by modulation of the coefficients of one or more regions of the said partition. 
[0088] Correlatively, the present invention also concerns a decoding device adapted to implement the decoding 
40 method according to the invention. 

[0089] The present invention also concerns a method of processing an image, comprising the following steps; 

inputting an image; 

inserting a first watermarking information into said inputted image; 
45 - applying a distortion to the image obtained in said insertion step; 

detecting said first watermarking information from the image obtained in said distortion applying step; and 
determining an amount of a second watermarking information which can be inserted into said inputted image in 
accordance with the result of the detection obtained in said detection step 

so [0090] It further comprises the step of displaying the amount of the watermarking information determined in said 
determining step. 

[0091] Correlatively, the present invention also concerns an image processing apparatus, comprising; 

means for inputting an image; 
55 . means for inserting a first watermarking information into said inputted image; 

means for applying a distortion to the image obtained by said insertion means; 

means for detecting said first watermarking information from the image obtained by said distortion applying means; 
and 
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means for determining an amount of a second watermarking information which can be inserted into said inputted 
image in accordance with the result of the detection obtained by sa.d detect.on means 

[0092] This apparatus further comprises the means for displaying the amount of the watermarking information de- 

l^S.^^l^^conc.rn. a method of inserting a watermarking signal which can take into account 
Se same tK^SLmined imperceptibility and de.ectabi.ity constraints in order to ensure invis.b.l.ty of the wa- 
termarking signal in the image and correct detection at the time of decoding. 

Iu0941 in the article entitled "A method for signature casting on digital .mages by I. PITAS, in Proc. iaP.W» 
2^5 2 8 SepTembe 1996, there is described a method of inserting a watermarking signa ,n a dig.ta I image ,n wh.ch 
fh weight cogent « is determined so as to ensure detectability of the signal with a fixed ***** 
00951 This weighting coefficient must be greater than a minimum value, which can be referred o as the detect.on 
L^lLe so as to a,low the detection of the inserted signal with the correct detection probability level. 
JSSHhE minimum va ,ue or detect.on amplitude depends notably on the length of the watermarking s.gna. and on 

" T^TwZTe^Z c^ent is fixed in the aforementioned document without taking info account the 
5, SS^^mL wh.ch guarantees invis.bil.ty of the watermark.ng signa, inserted >n the <^maga 
£Sl To th's end the present invention relates to a method of inserting a watermarking s.gnal .n a set of coeff c ents 
™ 6 n«Tnn a dStal image in which at least one subset of coefficients is modulated according toa linear model using 
^SSS^ing coefficient being greater than a min.mum value determined as a funct.on notab.y 

*ZZL£Z2Z ding to the .ength of the watermarking signa, ensuring imperceptibly of the watermarkmg ,s gna.. 
OlS Casting a maximum value for the weighting coefficient, which can be referred to as y.sual amp. ude^ 
makes it possible to take into account, when a watermarking signa. is inserted in a set of coeff , cents representing a 
diaital imaoe constraints related to the imperceptibility of the inserted watermark.ng signal. 
ra?01] Such an Son method makes ff possible to determine the conditions which make possible the insertion of 
L ^ lrmarktg signal in an image whilst meeting predetermined imperceptibility and detectab.l.ty criterion or on the 
oZ t££d2* -if. under certain conditions, the inserts of the watermarking signa. meets the predeterm.ned 

30 iXr™ d St^SE^ of imp.ementat.on of the present invention, the insertion method^ 
ndudes S ^determining compatible values of the length of the watermarking signal, of the we,ght,nE .coeff cient 
Lnd ErZJ ^robab Iffv evel so that the value of the weighting coefficient determined is less than or equal to said maximum 
value rndgTeate! than "r equal to said minimum value for said determined values of the length of the watermarking 

" Sfi" t?^^^^e^ .or the weighting coefficient, the conditions of detectability and 

thP watermarkino sianal and of the detection probability level in order to satisfy this inequality. 
.o Z^ ^2^Zp^erre6 characteristic of the invention, during the step of calculating the said maximum va.ue. 
a function independent of said set of coefficients to be watermarked is used. 

^^Ju-STmodel used for calculating the maximum value or the visual amplitude ,s independent of the image 

loi^TnT^terred embodiment of the invention, the coefficients are spatio-trequency "^^Jf^* 
4S LpatSfrequency transformation of this image, the ca.culat.on of said maximum va.ue of the weighting coefficient de- 

SSS? T^^Z^ZTwer of a watermarking signa, to app.y a spatio-f requency 
to°e image also referred to as transformation into sub-bands of the image, of the discrete cosine transformalon 
(OCT) ZZZ transformation type, so that the modulation for inserting the watermarking signa. ,s .mp.emented on 
so Qnatin-freauencv coefficients of a sub-band of the image. 

WMl Te inserton method according to the invent.on app.ies particularly we., to the conventional techniques of 
!nse^L by mo^uSon inserted in a transformed domain obtained by one o, the aforementioned spatio-frequency 

loiTor^practice, said calculation function used for calculating the maximum va.ue or visual amplitude a. as a 
ss function of the length P of the watermarking signal w is of the type: 
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25 



where 



a base (S) is a base value, dependent on the transformation used S and on the (base) sub-band considered for the 
insertion, of the maximum weighting coefficient ensuring imperceptibility when a single coefficient of said sub-band 
is modulated, 
10 p is strictly greater than 2, and 

E[lwlP] is the mathematical expectation of the function ItvlP. 

[0111] According to another preferred characteristic of the invention, said minimum value is determined by using a 
calculation function also depending on the variance of the coefficients to be modulated, on the distribution of the wa- 
rs termarking signal and possibly on the variance of an additive noise decorrelated from the coefficients modelling any 
distortions applied to coefficients. 

[0112] Thus the minimum value of the weighting coefficient or detection amplitude can be calculated by taking into 
account notably the distortions which the image may undergo, notably during normal processings consisting of com- 
pression of digital images, in so far as such distortions can be modelled in the form of an additive noise decorrelated 
20 from the coefficients representing the image. 

[0113] In practice, the said function of calculating the detection amplitude a D as a function of the length P of the 
watermarking signal w is of the type 



D JbP^c 



where a* is equal to the variance of the coefficients to be modulated, ajj is equal to the variance of the additive 
noise and a, b and c are constants which depend on the distribution of the watermarking signal and on the required 
30 detection probability level. 

[01 1 4] According to a preferred characteristic of the invention, the watermarking signal is a pseudo-random sequence 
with a predetermined distribution with a zero mean. 

[0115] The use of such a watermarking signal increases the robustness of the inserted watermarking signal with 
respect to an intentional attack aimed at extracting or modifying the watermarking signal, typically by statistical analysis, 
35 the kernel of the pseudo-random sequence produced from a predetermined distribution also being able to be given by 
a secret key. 

[0116] In a particularly advantageous embodiment of the invention, the method of inserting a watermarking signal 
includes the following steps: 

40 - calculating a maximum value of the weighting coefficients for a length of the watermarking signal equal to the 
cardinal number of the set of modulatable coefficients ensuring imperceptibility of the watermarking signal at the 
limit of perceptibility; 

calculating the length of the watermarking signal for a minimal value of the weighting coefficient equal to said 
maximum value calculated and for a predetermined detection probability level; and 
45 - comparing said calculated length with the cardinal number of the set of modulatable coefficients. 

[0117] Thus, by determining the maximum value of the weighting coefficient at the limit of perceptibility, that is to say 
in the eventuality of all the coefficients available being modulated, it is possible to determine the minimum length of 
the watermarking signal which makes it possible to obtain a correct detection probability with a weighting coefficient 
so equal to the maximum value calculated. 

[0118] The comparison of this length thus calculated with the cardinal number of all coefficients available in order to 
effect the insertion makes it possible to determine whether this insertion is actually possible with the required detection 
probability level. 

[0119] In addition, if the calculated length is less than the cardinal number of the set of coefficients, it can be certain 
ss that the imperceptibility criterion will indeed be met by modulating a subset of coefficients with a cardinal number equal 
to this calculated length, since the value of the weighting coefficient was calculated in the eventuality of all the available 
coefficients actually being modulated. 

[0120] In practical terms, if, at said comparison step, said calculated length is less than or equal to the cardinal 
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number of the set of modutatabie coefficients, a watermarking signal with a length at least equal to said calculated 
length is inserted by modulating a subset of coefficients with a cardinal number equal to said length, according to a 
linear model using a weighting coefficient equal to said calculated maximum value. 

[0121] Conversely, if at said comparison step said calculated length is greater than the cardinal number of the set 
s of modulatable coefficients, there is calculated the detection probability level obtained for a minimum value of the 
weighting coefficient equal to said maximum value calculated for a watermarking signal length equal to the cardinal 
number of all the modulatable coefficients. 

[0122] Thus the decision to insert a watermarking signal or not on the set of modulatable coefficients can be taken 
in the theoretical knowledge of the detection probability level obtained by default. 
10 [0123] In another practical embodiment of the invention, the insertion method includes a step of calculating a thresh- 
old value of the length of the watermarking signal determined so that said minimum value of the weighting coefficients 
is equal to said maximum value of the weighting coefficients for a length of the watermarking signal equal to said 
threshold value and a predetermined detection probability level. 

[0124] It is thus possible to determine the minimum length threshold which the watermarking signal must meet in 
is order to simultaneously meet the two criteria of perceptibility and detection when the watermarking signal is inserted 
in a digital image. 

[0125] Correlatively, the present invention also concerns a device for inserting a watermarking signal in a set of 
coefficients representing a digital image, in which at least one subset of coefficients is modulated according to a linear 
model using a weighting coefficient, the weighting coefficient being greater than a minimum value determined as a 

20 function notably of the length of the watermarking signal and of a detection probability level. 

[0126] This insertion device comprises means of calculating a maximum value of the weighting coefficient as a func- 
tion of the length of the watermarking signal ensuring imperceptibility of the watermarking signal. 
[0127] This insertion device has characteristics and advantages similar to those described previously since it is 
adapted to implement the insertion method according to the invention. 

25 [0128] The present invention also concerns a computer, a digital image processing apparatus, a digital printer, a 
digital photographic apparatus and a digital camera adapted to implement the method of determining a partitioning, 
the methods of inserting a watermark, the method of processing an image and the decoding method in accordance 
with the invention. 

[0129] These devices for determining a partitioning, inserting a watermark, processing an image and decoding, this 
30 computer, this digital image processing apparatus, this digital printer, this digital photographic apparatus and this digital 
camera have characteristics and advantages similar to those described with reference to the methods which they 
implement. 

[01 30] The present invention also relates to a computer program product and to an information storage means which 
can be read by a computer or a microprocessor, possibly removable, which stores a program implementing the method 
35 of determining a partitioning and/or inserting a watermark and/or processing an image and/or decoding in accordance 
with the invention. 

[0131] Other particularities and advantages of the invention will also emerge from the following description of the 
invention. 

[0132] In the accompanying drawings, given by way of non-limitative examples: 

40 

Figure 1 is an overall diagram illustrating a device for inserting a watermark into a digital image; 

Figure 2 schematically illustrates a step of spectral decomposition employed when inserting a watermark; 

Figure 3 is an algorithm of the method of determining a partitioning in accordance with a first embodiment of the 

invention; 

45 - Figure 4 is an example of partitioning of a region by division; 

Figure 5 is an algorithm of the method of determining a partitioning in accordance with a second embodiment of 
the invention; 

Figure 6 is an algorithm of the decoding method according to a first embodimenl of the invention; 
Figure 7 is an algorithm illustrating the application of the method of determining a partitioning of Figure 1 to a digital 
so image; 

Figure 8 is an algorithm illustrating the application of the decoding method of Figure 4 to a digital image; 
Figure 9 is a block diagram illustrating a variant embodiment of the method of determining a partitioning in accord- 
ance with a third embodiment of the invention; , 

Figure 10 illustrates the application of the method of determining a partition to a digital image according to a 
55 preferred embodiment of the invention; 

Figure 11 is an algorithm of the method of determining a partition in accordance with a first way ol the preferred 
embodiment of the invention; 

Figure 1 2 is an algorithm of the method of determining a partition in accordance with a second way of the preferred 
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embodiment of the invention ; 

Figure 1 3 is an algorithm of the insertion method according to one embodiment of the invention; 
Figure 14 illustrates, by means of a graph, the insertion method of the figure 13: and 

Figure 15 is a block diagram illustrating a device adapted to implement the methods in accordance with the inven- 
tion. 



[0133] A description will first of all be given, by reference to Figure 1 , of the insertion of a watermark into a set of 
coefficients representative of a digital image i. 

[01 34] This digital image I can be represented by a set of coefficients, either in the spatial domain or in a transformed 
10 domain, the coefficients in this latter case being hybrid coefficients, that is to say coefficients located both in the spatial 
domain and in the frequency domain. Such a representation of the image is obtained, for example, by using a sub- 
band originating Irom a space-frequency decomposition of the image, for example a discrete wavelet decomposition. 
[01 35] Here, the insertion of the watermark is achieved via a spread-spectrum insertion technique, by modulation of 
coefficients of a space-frequency representation of the image, obtained by the space-frequency transformation of the 
15 digital image I. 

[0136] By way of example, a discrete wavelet decomposition S can be used, represented diagrammatically in Figure 
2. This space-frequency decomposition is well known in the image-processing field, and only the principle thereof is 
reiterated below. It makes it possible to segment the image into sub-bands of frequencies and to obtain hybrid coeffi- 
cients, that is to say spectral coefficients which are also located in space, here in the plane of the image. 

20 [0137] The image I consists of a series of digital samples The image I is represented, for example, by a series of 
bytes, each byte value representing one pixel of the image I , which may be a black-and-white image with 256 grey levels. 
[0138] The multi-resolution spectral decomposition means consist of a sub-band decomposition circuit, or analysis 
circuit, formed by a set of analysis filters, associated respectively with dividers-by-two. This decomposition circuit filters 
the image signal I in two directions, into sub bands of high spatial frequencies and of low frequencies. The circuit 

25 includes several successive analysis blocks for breaking down the image I into sub-bands according to several levels 
of resolution. 

[01 39] By way of example, the image I is broken down here into sub-bands with a level of decomposition d equal to 3. 
[0140] A first analysis block receives the image signal I and filters it through two respectively low-pass and high-pass 
digital filters, in a first direction, for example horizontal. After passing through dividers-by-two, the resultant filtered 

30 signals are in their turn filtered by two respectively low-pass and high-pass filters, in a second direction, for example 
vertical. Each signal is again passed through a divider-by-two. Hence four sub-bands LL,, LH, , HL, and HH, are 
obtained at the output of this first analysis block, with the highest resolution in the decomposition. 
[0141] The sub-band LL, includes the low-frequency components in the two directions of the image signal I. The 
sub-band LH-, includes the low-frequency components along a first direction and high-frequency components along a 

35 second direction of the image signal I. The sub-band HL, includes the high-frequency components along the first 
direction and the low-frequency components along the second direction. Finally, the sub-band HK, includes the high- 
frequency components along the two directions. 

[0142] A second analysis block in its turn filters the sub-band of low frequencies LL, so as, in the same way, to supply 
four sub-bands LL 2 , LH 2 , HL 2 and HH 2 with intermediate level of resolution in the decomposition. A third analysis block 
40 then filters the sub-band of low frequencies LL 2 so as to supply four sub-bands LL 3 , LH 3 , HL 3 and HH 3 . Finally in this 
example, the sub-band LL 3 is in its turn analysed by a fourth analysis unit in order to supply four sub-bands LL 4 , LH 4 , 
HL 4 and HH 4 with the lowest resolution in this decomposition. 

[0143] Hence thirteen sub-bands and four levels of resolution are obtained. Obviously, the number of levels of res- 
olution, and consequently of sub-bands, may be chosen differently. 
45 [0144] Obviously, other types of spectral transformation could be used, such as the discrete Fourier transform, the 
discrete cosine transform or the Fourier-Mellin transform. 

[0145] In a general way, frequency sub-bands are obtained forming a set of spectral coefficients into which a water- 
mark can be inserted. 

[0146] The high-frequency sub-band of highest resolution HHj can be considered here in order to insert a watermark. 
so This sub-band HH-, thus supplies a set of coefficients X, equal to N in size, for example. 

[0147] Let X = {Xj, 1 < i < N} be the set of coefficients representative of the digital image I. 

[0148] Obviously, it would be possible to consider a subset of coefficients of this set X, representative only of a part 
of the digital image I. 

[0149] In order to insert a watermark into this set of coefficients X, one technique consists in inserting a pseudo- 
55 random signal, by spreading its spectrum : so as to render this signal undetectable by spectral or statistical analysis. 
[0150] Let w = (Wj, 1 £ j < P} be a watermark of size P < N, a pseudo-random signal with known distribution and 
average of zero. The most widespread distributions are the binary distribution {-1 , 1 }, the uniform distribution over [-1 , 
1] and a centred normalised Gaussian distribution N (0, 1). 
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[0151] The linear insertion formula is : 

XV = X- + b ctj w p with 1 < j < P, 

5 

in which {X jf 1 < j < P] is a subset of the set of coefficients X, b is an information bit, and is a weighting coefficient, 
also called modulation amplitude. 

[0152] The set {j, 1 < j < P} is also called information-bit insertion carrier. 

{0153] In order to insert a signal which may be formed from several information bits, it is worth partitioning the set 
io of coefficients so as to determine the number of regions or the number of insertion carriers available in this set. 

[0154] To do that, the insertion device 10 in accordance with the invention comprises a device for determining a 
partition 11 and insertion means 12 

[0155] The device for determining a partition 11 is suitable for carrying out adaptive partitioning of the set of coeffi- 
cients as a function of the detectability of an inserted watermark. 
is [0156] This method of adaptive partitioning consists in varying the number of coefficients modulated in order to insert 
an information bit of the watermark, so as simultaneously to satisfy the criteria which are set of invisibility and of correct 
detection probability. 

[0157] A description will first of all be given ; with reference to Figure 3, of a method of determining a partitioning into 
distinct regions of a digital image which makes it possible to maximise the number of regions obtained in order to 

20 maximise the number of information bits inserted in the image. 

[0158] Consideration is given to a representation of the digital image by a set of coefficients, either in the spatial 
domain, or in a transformed domain, the coefficients being in the latter case hybrid, that is to say located both in the 
spatial domain and in the frequency domain. Such a representation of the image is for example obtained by using a 
sub-band issuing from a spatio-f requency decomposition of the image, as a discrete wavelet decomposition previously 

25 described. 

[0159] The partitioning of the image is effected in order to insert a watermarking signal by modulation of the coeffi- 
cients of each region. 

[01 60] According to the invention, this partitioning into regions is performed by an adaptive partitioning as a function 
of a criterion of detectability of an information bit inserted on this region. 
30 [0161] The number of regions obtained under constraint of the detectability criterion makes it possible to define a 
maximum capacity of the digital image. This maximum capacity corresponds to the number of information bits which 
it is possible to insert in the digital image whilst complying with the detectability criterion. 

[0162] Initially, at a step E100, this capacity Q is initialised to 0 and the starting point is an initial partitioning which 
is limited here to a single region R of size P. 
35 [0163] Naturally, the initial partitioning could be different and notably include already several distinct regions with 
sizes which might be different. 

[0164] It is assumed that the size P of this initial region R is greater than a minimum size denoted P stat which corre- 
sponds to the minimum size of a statistically significant sample of coefficients for detecting a watermarking signal 
inserted in these coefficients. Typically, this minimum size can be fixed at 100. 
40 [0165] This initial region R can consist of the set of coefficients of the representation under consideration or a subset 
of this representation. 

[0166] A test step E101 then tests the detectability of the insertion of a watermarking signal in the initial region R. 
[0167] This detectability test can be effected by actual modulation and detection operations on the coefficients of 
the initial region R. 

45 [0168] In practice, all the coefficients of the initial region R are copied into a working memory so as not to directly 
watermark the image. 

[0169] Next all the coefficients of the initial region R are modulated by inserting the watermarking signal. 
[0170] Conventionally, the following procedure is followed in order to effect this modulation: 
[0171] The initial region R consists of a set of coefficients: 
so X: {Xj, U 1. ...P) 

[0172] Normally a linear modulation is used of the type: 

X* = Xj ± ttjWj, with i= 1 P 



55 



where W : {W j: i = 1 , ... : P) is the watermarking signal generally chosen as a pseudo-random sequence with a 
known distribution and a nil mean. 

[0173] The value of the sign ± of the modulation depends on the binary value to be inserted: for example the sign - 
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corresponds to the value 0 and the sign + corresponds to the value 1 

[0174] It is possible to choose, by way of example, a pseudo-random sequence W which follows a uniform law on 
the interval [-1, 1]. 

[0175] Naturally, any pseudo-random watermarking with known distribution and nil mean can be suitable. The most 
5 usual distributions for the watermarking signal W are, apart from the aforementioned uniform distribution [-1, 1]. the 
binary distributionM.1^ and the centred standardised Gaussian distribution N(0, 1). 

[0176] The modulation can optionally be protected by means of a secret key characterising for example the kernel 
for reproducing the pseudo-random sequence W. 

[0177] The term ctj is the modulation amplitude, or weighting coefficient, applied to the coefficient Xj. 
io [0178] A constant weighting coefficient can be used for all the coefficients of the region R, such that <xj = a v for any 
i, the value of the weighting coefficients being equal to a maximum value ensuring imperceptibility of the watermark- 
ing signal for the modulation of the set of coefficients representing said digital image. 

[0179] This maximum value a v (referred to as Just Noticeable Difference or JND amplitude) decreases in fact with 
the length of the sequence of the watermarking signal W, and consequently with the number of modulatable coefficients. 
is [0180] Given that the capacity Q will be maximised when all the coefficients of the representation ol the digital image 
are used for the modulation, it is judicious to determine a modulation JND amplitude for a number of modulatable 
coefficients equal to the cardinal number of the set of coefficients representing the image. 

[0181] This maximum value a v or JND amplitude corresponds to the maximum modulation amplitude which can be 
used in the aforementioned linear insertion model beyond which an observer is capable of visually detecting a change 

20 in the reconstructed watermarked image. 

[0182] In order to determine this JND amplitude, it is possible to use a visibility model which makes it possible to 
predict the visibility of a watermarking operation according to different parameters, namely the representation of the 
signal through a spatio-frequency transformation S used, the sub-band considered for the insertion, the type of distri- 
bution of the sequence of the watermarking signal W and the length P of the sequence W. 

25 [0183] A simple model, developed by WATSON and described in the article entitled 'Visibility of wavelet quantization 
noise", A.B. WATSON et al, IEEE Trans, on Image Process, 6(8), 1164-1175, 1997, makes it possible to predict the 
visibility of a set of modulated coefficients from the measurement of visibility of a single modulated coefficient. Reference 
can advantageously be made to this document for the detailed description of this model. 

[0184] Thus use is made of a function dependent on the length P of the watermarking signal W, the type of trans- 
30 formation S used and the (base) sub-band considered for the insertion, but independent of all the coefficients to be 
watermarked X. 

[0185] This calculation function used can be written: 

OL has ~(S) 

p p (Eii^nr 

where 

40 ct t>ase( s ) ' s a Dase value, dependent on the transformation used S and on the (base) sub-band considered for the 

insertion, of the maximum weighting coefficient ensuring imperceptibility during modulation of a single coefficient 
of this sub-band, 
|i is strictly greater than 2, and 

EllwlP] is the mathematical expectation of the function livlP. 

45 

[0186] The base values a base (S) can be measured once and for all, for each sub-band of coefficients in the wavelet 
decomposition, from a single psychovisual measurement, and be stored in a table of visual amplitudes. 
[0187] (5 is the exponent of a Minkowsky sum and can be chosen equal to 5, for example. 
[0188] The mathematical expectation £[lwlP] corresponds to an estimation of the mean of the function livfP. 
so [0189] This visibility model does not take into account the image I to be watermarked itself and is independent thereof. 
This is equivalent to considering that the image I is unilorm. It is a question of a "worst case" model since the presence 
of the image signal makes it possible to visually mask the modulation itself. 

[0190] The JND amplitude can thus be calculated for a length P of the watermarking signal W equal to the cardinal 
number of all the modulatable coefficients of the sub-band under consideration. 
55 [0191] However, it is judicious to exploit the fact that the image signal itself makes it possible to mask the modulation. 
[0192] Thus, in this example, use is made, for each coefficient X ( to be modulated, of a weighting coefficient in 
accordance with a law of the form = kj.ay, where kj is a modulation factor dependent on the coefficients situated 
close to the coefficient considered on the region R and cx v is equal to the maximum JND value of the weighting coef- 
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IS Thus each coefficient wi.. be modulated according to the local content, which maKes it possible ,o locally 
Increase the amplitude of the modulation , to « * L'deteSbi.ity amplitude is next ca.cu.ated from the de- 

mined threshold value denoted *■ , ft coefficients X' and me mark*!, signal W is 

SKL'iSSr. 2 SEES. 2™a,; dean watermartting Che .mage vrft.n the reee,, - ~ 
S ,:Sl m itSL opr. W, k gr.ate, ,han the P^«-™^ ^J, „,. one desended in ,h. anicle 
10196] A Mandated statistical last can also ba used for s ™ ? ° 215 . 218 carter ,996. 

spending to a fixed detection probability level, for «aj^j»£ 5 *. weighting coefficient a, it 

detectability criterion on this region R. . nenHinn on the sjze P ot lhe region R, the detection 

S5=S=SSra«ri=is==: 

function: 

Pmin,X) = (|e) ! |dd%.cand| 

.^^^ 

S°Sn SingCaaf o, coefficients X demg modulated. Ihis tunatxm also depende on Una vananc. o, 
S"r n m aSr.^ ~™,o anovide ,o, an y distortions »ni* Ida .mage ma, undergo ana, mo*,—, dee 
tor example to a compression of the image. ^« rr oi a t«ri from the coefficients it is possible to use 

modulation. onc.irp invisibilitv of the inserted watermarking 

coefficients X there also corresponds a minimum modulation length P min (X) wnicn manes w 
;S^ M " y o°i S„"Srn'da,disad hypothes* lasl. m,a minimum .englt, P^X, re man det.neo d, 
the function: 

P min(X) = (SZf j* W*\ * « 8 <4 ♦ o*> ♦ ca 2 ♦ d] 
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50 where a, b, c and d can be assimilated to constants. decoding tochoose a thresh- 
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of the digital image without having to transmit to the decoder the partitioning effected at the time of insertion of the 
watermarking signal W. 

[0208] In this example embodiment, a partitioning by iterative division of the digital image is effected. It is conse- 
quently considered, at step E101. that the detectability criterion is satisfied when the calculated detection amplitude is 

s greater than or equal to the predetermined threshold Tc. 

[0209] If the detectability criterion is not satisfied, then it is probable that the detectability criterion is also not satisfied 
in subregions of the initial region R of smaller size because the detection probability generally decreases with the size 
of the regions. A step E102 of stopping the method of determining a partitioning makes it possible to stop the process 
of iterative division of the image. 

io [0210] Insertion not being judged possible in the initial region R, the capacity Q of the image is therefore nil. 

[0211] If on the other hand the detectability criterion is satisfied, that is to say here that the detection amplitude is 
greater than or equal to the threshold value Tc, it is possible that it is also satisfied for subregions of the initial region R. 
[0212] It is therefore judicious to attempt to divide the region R as long as the detectability criterion is satisfied in 
order to increase the capacity Q. 

75 [0213] There are of course several possible means for effecting the spatial division of the region R into subregions 
Rj, of sizes Pj, with j = 1, K. 

[0214] In this example a division into a quadtree is chosen, as illustrated in Figure 4. 
[0215] This division consists of dividing each region into four subregions of identical size. 

[0216] Thus, by way of example, in Figure 2, a digital image of initial size N is being partitioned into two sub-regions 
20 R5 and R6 of sizes N/4, seven subregions R1 , R2, R3, R4, R7, R8, R9 of sizes N/16 and four subregions R10, R11 , 
R12. R13 of sizes N/64. 

[0217] Thus, at the division step E103, the region R is divided into subregions Rj, j = 1 K of sizes Pj. In a case of 
a division into a quadtree, K = 4 and Pj = P/4. 

[0218] At step E104 a first subregion of index j (for example j = 1) is considered and a value of potential gain in 
2S capacity denoted G Q is initialised to - 1 . 

[0219] The partitioning determination method then includes a prior step E 105 of comparing the size Pj of the region 
Rj with the minimum size P stat corresponding to the minimum size of a statistically significant sample for the detection 
of a watermarking signal inserted in said region as described above. 

[0220] If the size PjOf the region Rj is less than this minimum size P slat , at step E108the index J = J+1 is incremented 
30 in order to process the following sub-region Rj +1 . 

[0221] Otherwise, the detectability criterion is checked on the subregion Rj in the same way as previously described 
with reference to step El 01. 

[0222] In practice, the watermarking signal is inserted on the coefficients of the said region Rj, a detection amplitude 
is calculated and this detectability amplitude is compared with the threshold value Tc. 
35 [0223] The detection amplitude can be a correlation value between the modulation signal X* and the carrier of the 
watermarking signal W or the result of a standardised hypothesis test. 

[0224] If the detectability criterion is verified, the value of potential gain in capacity Gq is increased by one unit at 
step E107. 

[0225] Then at step E108 the index j of the subregions is incremented and it is checked at the test step E109 that 
•*o this index j is less than or equal to the number of K. 

[0226] In the affirmative, steps E105 to E109 are reiterated on the following subregion. 

[0227] Otherwise, when the index j = K+1, it is checked in a test step E110, whether the potential gain in capacity 
value Gq is strictly positive, that is to say whether there exist at least two subregions Rj of the initial region R for which 
the detectability criterion is validated. 
•fs [0228] In the negative, the partitioning of higher rank carried out on the initial region R is not validated and this region 
R is preserved. Such is the case when the detectability criterion has not been validated for any of the subregions Rj 
or only for a single subregion Rj. 

[0229] In the latter case, it is preferable to keep the region R for inserting the watermarking signal since it is ol size 
P greater than the size of the single validated subregion Rj. 
so [0230] At step E112, the value of the capacity Q is then increased by one unit and the region R of size P in the 
partitioning is kept. 

[0231] An actual insertion step E113 by modulation of the coefficients of the region R can then be implemented in a 
normal fashion using the linear insertion model described above. 

[0232] An elimination step E114 then makes it possible to eliminate the region R from a stack of regions to be proe- 
ms essed. 

[0233] On the other hand, if at step E 1 1 0 the value of the potential gain capacity Gc is strictly positive, in an addition 
step E111 all the subregions Rj for which the detectability criterion is validated in the stack of regions R to be processed 
are added. 
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10234] A return step E115 makes it possib.e to reiterate a,, the steps E103 to E113 in order to process the regions 
of the stack in a recursive fashion t , e „ ri . H a ran a r itv value Qm corresponding to the 

,„, coefficients d .he di«e,en. regions gained by .he P'» to "^e m e^r.«»v 

E113 identical to those of the method of determining a partitioning of the embod.ment of F,gure 3 and wh.cn 
^Tas^ 

20 [^r^he method includes a comparison step E116 in which this value of the capacity Q is compared with the 

l^ZZTZTJ^ a panning a.so includes, in this ^T^lT^T^ ^ 
S regions R to be processed in the stack by detectabi.ity amplitude value calculated on these eg»n _R_ 
30 [0246? Thus a return step E11S reiterates the partitioning from the division step E103 for a reg.on R .n the stack 

having the highest delegability amplitude amongst all the other regions { ° be J> rooessed .. . soastODroceS s 
raw? Asortingvariantcanconsistofsortingtheregionstobeprocessed.ndecreas.ngorderofs.zesoastoprccess 

as a priority, the largest regions which a priori provide better delegability. 

[0248] Naturally, a partitioning by iterative merging of the d.g.tal image could be earned out in the same way, 
55* £ S^gC^o^ with smaller acceotab.e sizes, greater than or e.ua, to the 
K U ^For e eac e h d ^^p2lSc^ tor which the de.ectabi.ity criterion is not validated, this region is merged 
£CU is if possible a region for wh.ch the detectabi.ity criterion is no, va.idated 
^ a 5 S 2] t0 SSS region can be a region satisfying the detectabi.ity criterion so as to make detection 

[SSt There are many possibi.ities for the choice of the reg.on to be merged. It is f^^^SZ 
Son the detection amplitude calculated on each of the regions and merge for example the reg.ons associated with 

loasT?^^ 

and consists of deciding whether or not fou. adjacent blocks, for example R1. R2. R3 and R4 ,n Figure 4, mus 

so A description will now be given, with reference to Figure 6, of the ^ 

t possib.e to find the partitioning when the watermarking signal inserted ,n a d.grtal "^VT^^ fae known a 
[0256] indeed, because of the adaptation of the partitioning to the d.g.tal .mage, this partrt.on.ng w.ll not be known a 

jS^]" The' deling method then consists, in genera, terms, o, testing a., the acceptable ■^JjJ^ 
ss Iterative imp.ementation technique is used as the one used when the partit.on.ng • determ.ned for .nsert.ng .he wa 

;oS8rTura a desch P tion II be given of an example embodiment in the case of a segmentation by dMdhgby 
££ in a "ua^re 7ymme,nca. to'tha, using the partitioning determination method of the embod.ment described 
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with reference to Figure 3. 

[0259] At an initialisation step E200, there is an initial partitioning limited to a single region R of si7e P. This region 
consists of all the coefficients of the same representation of the digital image as the one used at the time of insertion 
of the watermarking signal. 

5 [0260] A test step E201 makes it possible to verify that the size P of the initial region R is indeed greater than the 
minimum size P slat as described previously. 

[0261] In the affirmative, at step E203 a criterion of detectability of the information bit inserted on the region R, 
obtained by demodulation of the coefficients on this region R, is verified. 

[0262] In practice, a detection amplitude is calculated for all the coefficients of the region R and the detectability 
10 criterion is validated by comparing this detection amplitude with a predetermined threshold value Td for decoding. 
[0263] It is possible to proceed by means of a standardised statistical test, and the result compared with the prede- 
termined threshold value Td. 

[0264] As explained before, it is necessary, if the same detectability test is used during coding and decoding, to 
choose a threshold value Tc on coding greater than that Td used at decoding if it is wished to be able to find the same 
15 partitioning on decoding of the image even when the latter has been made noisy by post-processing. 

[0265] This condition thus makes it possible to find the same detectability criterion which constrains the partitioning 
of the digital image without having to transmit, to the decoder, the partitioning effected when the watermarking signal 
W was inserted. 

[0266] If the detectability criterion is not verified at step E203, a division is carried out, at step E204, of the regbn R 
20 into subregions Rj and sizes P| with j = 1 , .... K. Then the same partitioning structure is used as the one used on coding, 
that is to say, in this example, a division of the region R into four subregions R, of identical sizes. 
[0267] Then, in an addition step E205, the subregions Rj are added in a stack of regions R to be processed and, at 
a return step E202, the decoding method is reiterated recursively as from the test step E201 on each region R of the 
stack. 

25 [0268] When at the step of calculating the detectability criterion E203, the latter is verified for a region R, there is 
extracted at the extraction step E206 the information bit inserted in this region R as an element of the method inserted 
in the digital image I. This extraction is effected in a conventional fashion by calculating a detection amplitude. The 
result of the test gives, in absolute value, a detection amplitude to be compared with the threshold value Td of the 
decoding, the sign giving the value 0 or 1 of the information bit inserted. 

30, [0269] It should be noted that this extraction is already partly effected during the detection step E203 so that the 
detection E203 and extraction E206 steps are structurally coupled. 

[0270] An elimination step E207 is implemented in order to eliminate the region R from the stack of regions to be 
processed, which was formed at the addition step E205, and the return step E202 then makes it possible to implement 
the decoding method recursively on each region R of the stack to be processed. 
35 [0271] A description will now be given, with reference to Figures 7 and 8, of a practical application of the method of 
determining a partitioning of a digital image in a particular embodiment of the invention. 

[0272] In this application, the watermarking signal is inserted by means of a technique of insertion by spectrum 
spreading, by modulation of coefficients of a space-frequency representation of the image, obtained by a spatio-fre- 
quency transformation S of the digital image I as previously described with reference to figure 2. 
40 [0273] In general terms Sj sub-bands are obtained, with J = 1 , M and M equal here to 1 3 at a transformation step 
E300. 

[0274] Next, for each sub-band Sj of index j considered at step E301, a partitioning is effected at step E302 by 
implementing the method of determining a partitioning described previously to with reference to Figure 3. 
[0275] In this way a capacity Qj is derived tor each sub-band Sj. 
45 [0276] A step E303 of insertion by modulation of the coefficients next makes it possible to insert Q] information bits 
in the Qj regions validated. 

[0277] In steps E304 and 305, the following sub-bands are next considered as long as all the sub-bands Sj have not 
been processed, and the steps E302 and E303 of partitioning and inserting a watermarking signal of capacity Qj are 
reiterated on each sub-band Sj. 

so [0278] Next, an inverse transformation S' 1 is applied to the image in order to obtain the watermarked image, of total 
capacity Q equal to the sum of the capacities Qj calculated for each of the sub-bands Sj. 

[0279] Thus, each sub-band of the representation is considered independently for the adaptive partitioning which is 
carried out both at coding and at decoding. 

[0280] As illustrated in Figure 8, the decoding method in this embodiment is implemented in a similar fashion using 
55 the decoding method described with reference to Figure 6, the initial region R here corresponding each time to a sub- 
band of the decomposition of the image. 

[0281] A step E400, identical to step E300 : divides the image I into Sj subbands with j = 1 , ... M. 

[0282] A first sub-band Sj is considered at step E401 . The order of processing of the sub-bands is identical to that 
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used durinq coding in order to extract the inserted message in order. 

[02831 The adaptive partitioning is implemented at the partitioning step E402 in order to calculate the capacity q of 
the sub-band S, Next, at the extraction step E403. the Q, information bits are extracted by demodulat.ng the coefficients 
of the Q; reqions of the partitioning. . 
Eo2841 At the test steps E404 and E405, it is checked in the usual fashion that all the sub-bands have been processed 
and, in the negative, the decoding steps E402 and E403 (partitioning and extraction) are reiterated on the rema.n.ng 

[WW] n Tn this way, at the output ot the decoder, a message of Q bits is obtained, Q being equal to the sum of the 

capacities Q. of each sub-band Sj. ^~~-h 

[0286] In a variant of this method illustrated in Figure 9, it is possible to effect a partmomng on a post-processed 
image signal in order to obtain an a priori robustness of the watermarking signal to any postprocessing undergone by 

ISsTl^The principle of this variant consists of performing a delegability test in a test step E309 which is no longer 
coupled to the steps of partitioning E302 and insertion E303 described previously orWm£jH 
[0288] After the insertion of the watermarking signal in a region, a reverse transformation is first of all performed 
at step E306 in order to find the watermarked image. Then, in a step E307, post-processing is applied, such as a 
comDression of the digital image. . ... 

r02891 Next in a step E308, the transformation S is reiterated in order to find the watermarked coefficients, possibly 
made noisy by the post-processing carried out on the image at step E307, and the delegability criterion ,s calculated 

20 as previously described. , . 

r02901 Steps E300 lo E309 are reiterated on another acceptable partitioning of the .mage, here for example by 
dividing the relevant region R into four, and the detectability test is performed at step E309 on the four sub-regions^ 
[0291] This partitioning is adopted only if it represents a gain in capacity, that is to say, in this example of quadtree 
division if the detectability criterion is validated on at least two sub-regions. 

25 [0292] This variant embodiment makes it possible to directly take into account the post-proccss.ngs wh.ch w.ll be 

aDDlied to the digital image. ... . . 

[0293] A method of determining a partition according to a preferred embodiment of the mvent.on will now be de- 

scribed, by reference to Figure 10. 

r02941 The method of inserting a watermark is employed on a digital image I 

[0295] This insertion method first of all includes a method of determining a partition of a set of coefficients represent- 

[0296]' ' Thfs' plrtSs 5 applied as a function of a criterion of detectability of an information bit inserted on each region 

f0297] Pa It should be noted that, in practice, a partition of a digital image into a maximum number of regions satisfying 
a detectability criterion can be sought so as to ascertain the practical capacity of this image with regard to he ^sertion 
of a watermark. This maximum capacity of the digital image corresponds to the number of information bits wh.ch it is 
possible to insert into the digital image while complying with the detectability criterion. 

[0298] It is also possible, for a given watermark, comprising a given number of information bits of less than the 
maximum capacity of the image I, to seek the best possible partition of the image I which makes it P^le to insert 
this watermark of given length, for a given detection probability, and in regions of the largest possible sizes so as best 
to ensure invisibility of the inserted watermark. 

[0299] The method of determining a partition first of all includes a step of spectral decomposition E500 of the wavelet- 
decomposition type described previously (Discrete Wavelet Transform or DWT) 

F0300] Thus a set of spectral coefficients X is obtained, representative of the dig.tal image I. These spectral coeffi- 
cients are distributed into frequency sub-bands as illustrated diagrammatically in Figure 2, on completion of the de- 
composition step E500. The size of the set of coefficients X is equal to N = M x M for a square image. Obviously, the 
same method could be applied to a rectangular image. ^..A^Hnm 
[0301] In this embodiment of the invention, the melhod includes a step E501 of generating a centred pseudo-random 
sequence of size N equal to the cardinal number N of the set of coefficients X representative of the image I. 
so [0302] Let w = {w mn , 0 < m, n < M} be this pseudo-random sequence. 

r03031 This pseudo-random sequence w is formed from centred pseudo-random sub-sequences w*. 
[0304] The pseudo-random sequence can be seen as the merging of the pseudo-random sub-sequences w*. 
[0305] Here, the centred pseudo-random sequence w is formed from pseudo-random sub-sequences w* of identical 
size which correspond to a block, called unit block, of the set of coefficients X 
ss ro3061 Obviously the pseudo-random sub-sequences w k could differ in size from one another. 

0307] Here the size of these unit blocks corresponds to the minimum number of coefficients suitable for forming a 
region of a size which is statistically significant for satisfying the detectability criterion over this region 
[0308] As described later, when a standardised statistical detection test is applied to a set of coefficients, a m.n.mum 
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number of coefficients is necessary for this test to be significant. 

[0309] The pseudo-random sub-sequences w k are of a size A = L x L which is greater than or equal to this minimum 

number of coefficients. ... 

[0310] For preference pseudo-random sub-sequences w k are chosen to be equal in size to this minimum size. This 
5 characteristic is particularly advantageous when partition is achieved by iterative combination by combining blocks 

from an initial segmentation of the image into unit blocks of size L x L. The initial segmentation then corresponds to 

the maximum partition of the image which can be applied, beyond which the blocks of coefficients generated are no 

longer statistically significant for the insertion and detection of an information bit. 

[0311] For example, L = 8 is chosen. 
10 [031 2] The pseudo-random sequence w is then the merging of N/64 pseudo-random sub-sequences w k forming unit 

blocks 

[031 3] On each unit block, in accordance with the invention, a centred pseudo-random signal of predefined distribu- 
tion (Gaussian uniform, etc.) has to be generated on the basis of a function depending on a secret key K and on the 
index k of each unit block. This key can be denoted f(K, k). It is then necessary to ensure that, in practice, the average 
is of the signal w k is exactly equal to zero over each unit block. 

[0314] A tirst method for generating such a centred pseudo-random sequence consists in centring the sequence in 
a deterministic way, by symmetrisation. This method is suitable for generating a pseudo-random sub-sequence of even 

[0315] In this case, each centred pseudo-random sub-sequence w k is created by generating one half, A/2, of the 
20 pseudo-random samples of the sub-sequence via a pseudo-random number generator of known centred distribution, 
and by symmetrising the pseudo-random samples thus generated in order to obtain the other half A/2 of the samples 
of the sub-sequence w k . 

[0316] For a single-dimensional sequence, a set of samples of size A/2 {w jt 1 < j < A/2} is therefore generated, and 
this set is made symmetric in order to obtain the other half of the samples {Wj = -w A .j, A/2 < j < A}. 
25 [0317] In order to reinforce the random character of the sequence thus obtained, permutation can be carried out on 
the set of samples obtained {w^ 1 < j < A} by using a secret key. 

[031 8] A second method of generating a centred pseudo-random sequence consists in creating each centred pseudo- 
random sub-sequence w k by generating pseudo-random samples via a pseudo-random number generator and by 
redistributing the sum of these pseudo-random samples generated over each sample. 
30 [0319] Thus, first of all, by using the secret key f(K, k). a sequence of pseudo-random samples of the same size A 
as the unit block, i.e. {Wj, 1 < i < A}, is generated. 
[0320] Nexl the exact sum of these samples is calculated: 



35 



[0321] Then this sum S k is redistributed over the samples so as to obtain a centred pseudo-random sequence, that 
40 is to say one with a sum exactly equal to zero. 

[0322] Thus w* = (Wj-S^A, 1 < i < A} is obtained. 

[0323] The pseudo-random sequence w is thus composed of the merging of the centred sub-sequences w* in such 
a way that the partitioning of a set of coefficients modulated by the centred pseudo-random sequence w gives rise to 
subsets of coefficients also modulated by centred pseudo-random sub-sequences. 
45 [0324] Obviously, other techniques could be used to generate centred pseudo-random sequences. 

[0325] As illustrated in Figure 1 0, a step E502 of modulation of the set of coefficients by the centred pseudo-random 
sequence w is then implemented so as to insert an information bit onto the set of coefficients originating from the 
speclral decomposition. 

[0326] In practice, the set of coefficients is copied into a working memory, for example a random-access memory in 
so a computer, so as not to mark the image I directly. 

[0327] The same information bit is modulated onto the set of coefficients, for example b = 1 here, according to the 

modulation formula: 



55 



X'j = Xj + a- Wj, withO < i < N. 

[0328] The weighting amplitude a x is chosen in the usual way so as to guarantee invisibility of the inserted information 
bit. 
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the set of coefficients representative ol the digital image . jt DOSsible t0 maS k the modulation 

the watermark for the modulation ol the set of coefficients makes ft sjb , e 

■ [0332] Thus, eachcoefficient will be modulated as a funct.on o he ^J^SSv 

locally to augment the amplitude of the modufator panning to a possessed 

- L - P — Und6r90ne 
by the image I. eery} tn the set of modulated coefficients. 

tortion to be applied to each level of the partition. anolvina a distortion includes the following 

[0336] This is all the more advantageous since, here, the step E503 of applying 
20 sub-steps as illustrated in Figure 10: 

. inverse spectra, transformat.on E503a of the set of coefficients K, modulated in order to obtain a watermarked 

- ^cation proper, E503b, of distortion to the wat e ™ arked J* ^ a set Q , mooulated coefficients after 
. spectral transformation E503c of the watermarked image in order to obta.n a oi moo 

distortion. 

40 [0341] The means for displaying to the users Q bit ^ ?££2y descried. 

0342 Various techniques for iterative f 9™ ntat ™ ^ by reference to 

0343] A first way of performing the ^ n J^°^lt ^^tasis of an initial segmentation of the set of 
. Figure 11, In which the partitioning is achieved by ^^T^^^^a^ over these regions, 
coefficients into regions of a statistically sign, ican \ s,ze ^ a '^ 1 f^g^ modulation E502 and distortion E503 
[0344] Fo.lowing on from the steps of transformation E ™J™™™^ 5 £*Zl sel o1 coefficients to be parti- 
es described above, a first frequency sub-band will be considered, in a step E504, 

iSs, The initial segmentation o, this set of coefficients corresponds to a spatial segmentate o. the image I into 

ttJZZttTZXZZ. iteratLy, giving the notion o, 

SSTSS^ parthion takes place frequency sub-band by frequency sub-band of the signal arising from the 
[0 P iT SEES, is considered here, in a step E505, in an order of scanning o, the sub*and, for example from 
IS" 9 A :2£2T.5L over a firs, unH block wh.h has no, ye, been processed 
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[0351] In a calculation step E506, a statistical detection test is applied to this unit block of modulated coefficients XV 
[0352] An example ot a standardised statistical test is given below. 

[0353] Let 1 < i < P be the set of coefficients to be tested after insertion of the pseudo-random sequence w and 
application of p re-distort ion. 

[0354] The test consists in calculating, for this set of size P: 

Jv 

in which M represents an average value of correlation between the values of each coefficient X'j and of modulation 
Wj, and V the variance of these values. 
[0355] Hence in practice: 



P i = i 



20 and 

25 r —\ J = 1 



[0356] The calculation of this test value T requires knowledge of the pseudo-random modulation sequence w which 
is generated once and for all at the above-described generation step E501 
30 [0357] The detectability is given by the comparison of this test value T with a certain threshold P p , which corresponds 
theoretically to a probability of correct detection of p% if the test T follows a Gaussian law with variance equal to 1 for 
the two hypotheses tested (presence or absence of watermark). 

[0358] A comparison step E507 makes it possible to compare this test value T with the threshold value P p or correct 
detection rate. 

35 [0359] On completion of this comparison, if the calculated value T is higher than the correct detection rate, the block 
tested is held in memory, in a memory-storage step E509, as being a region available for inserting an information bit. 
[0360] In practice, a capacity-calculating step E51 0 can be implemented in such a way that, for every block of coef- 
ficients tested which is suitable for having an information bit inserted, the value of the capacity is increased by one unit. 
[0361] In contrast, if the calculated value T is less than the correct detection rate P p , a block of higher level is selected, 

40 in a combining step E508, that is to say, here : a parent block corresponding to the combination of four unit blocks, and 
the steps E506 ol calculating the test value T and of comparison E507 with the correct detection rate P p are reiterated 
over this new block. 

[0362] During this combining of child blocks in order to form a parent block, those child blocks are preferably chosen 
tor which the detectability criterion is not satisfied. 
45 [0363] For application of this detectability criterion, any standardised statistical test can be used for the detection, 
such as that described in the article entitled "A method for signature casting on digital images" by I. PITAS, in Proc. 
ICIP, pages 215-218, September 1996, which makes it possible to characterise the detection in terms of probability 
[0364] A threshold value P p can be chosen, corresponding to a fixed level of detection probability, 99.95% for ex- 
ample. 

50 [0365] Hence, for a given region defining a set of coefficients X of size P and with a fixed weighting coefficient a, it 
is possible to calculate a detectability amplitude and to compare the latter with a threshold value P p so as to validate 
or not validate the detectability criterion over this region. 

[0366] In a general way, the detection probability decreases with the size of the regions tested, such that, by com- 
bining regions in order to increase the size of the block ol coefficients used for inserting a watermark, the value of the 
55 test T is increased over this new block. 

[0367] Hence, on completion of these steps of adaptive partitioning E504 to E508, an image insertion capacity can 
be determined, equal to O bits, as can an appropriate partition of the image corresponding to this capacity Q. 
[0368] A watermark can then be formed, with length equal to O bits, and can then be inserted, bit by bit. into the 
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various reqions of the partition. . 
T03691 Obviously the step E520 of insertion proper of the watermark can be coupled w.th the determination of the 
partition in such a way that each region available on conclusion of the detectability test E506. E507. can actually be 
watermarked by the insertion of an information bit by modulation of the coefficients of this region. 
I-03701 The insertion step E520 makes it possible actually to insert an information bit of the watermark, by choosing 
a coefficient value b = 1 or = -1 in the above-described modulation formula, the alternating of the + /-1 bits conveying 
the content of the message inserted 

[0371] A centred pseudo-random sequence of known distribution can be used, different tor each reg.on of the par- 
tition, and initialised by a secret key K'. as illustrated in Figure 1 

[0372] It should be noted that this secret key K' may, if appropriate, be identical to the secret key K used dur.ng the 
qeneration of the pseudo-random sequence w for determining the partition. 

[0373] Moreover, a centred pseudo-random sequence w identical to that generated at the generation step E501 

f0374] be TpTrtTtion of an image I is shown diagrammatically by way of example in Figure 10, on conclusion of steps 
E504 to E508, in which, for example, 25 separate regions are obtained making it possible to insert a 25-brt watermark. 
[03751 Obviously, any other partitioning technique could be used. 

[0376] A second way of performing the preferred embodiment of the invention will now be described, by reference 
to Figure 12, in which the partitioning is achieved by iterative division of the set of coefficients 
103771 On conclusion of the steps of transformation E500, generation E501 , modulation E502 and distortion E503 
as described above, a first frequency sub-band will be considered, in a step E51 1 , as set of coefficients to be partitioned^ 
[0378] A quadtree-mode division will be considered here, in which one block, called "parent block" is always divided 
into four blocks of equal size, called "child blocks". 
[0379] Obviously, partition takes place frequency sub-band by frequency sub-band of signal originating from the 

spectral decomposition. 

25 [0380] Hence, in a step E51 2, an initial block corresponding to a first frequency sub-band is considered. 
[0381] This initial block is stored in memory in a list L of blocks to be processed 

[0382] Next, in a step E513, a first block B to be processed, coming from the list of blocks to be processed L, is 
considered. 

103831 A detectability criterion is applied on this block B in a checking step E51 4 

[0384] In practice as previously described, the value T of a statistical detection test on this block of modulated 
coefficients is calculated in a calculation step E51 4a. This step is analogous to the calculation step E506 described by 

103 e 8S] nC Vhe F detectability is then given by the comparison of this test value T with a certain threshold P p . which theo- 
retically corresponds to a correct detection probability of p% if the test T follows a Gaussian law with variance equal 
3S to 1 for the two hypotheses tested (presence or absence of watermark). 

103861 A comparison step E51 4b makes it possible to compare this test value T with the threshold value P p or correct 
detection rate. This comparison step E51 4b is analogous to the comparison step E507 described by reference to Figure 

[0387] On completion of this comparison, il the calculated value T is less than the correct detection rate P p , a check 
40 is made, in a test step E515, as to whether there are blocks remaining to be processed in the list L. 
[0388] If not, the recursive partitioning algorithm is terminated. 

[0389] Otherwise the following block B of the list L is considered at step E513, and the checking step E514 is reit- 
e rated 

[0390] It on completion of Ihe comparison step E514b, the calculated value of the test T is higher than the correct 
4S detection rate P D . detectability is tested over each child block of this parent block B in a checking step E51 6 

[0391] m practice, the parent block B is divided into four child blocks and the detectability criterion is checked on 

each of these blocks according to a step analogous to the previously described checking step E51 4. 

[0392] A test E517 makes it possible to verify whether the detectability criterion is valid over at least one child block. 

[0393] If yes. a division step E518 is implemented in order to actually divide the block B into four child blocks, and 
so an addition step E519 allows these child blocks to be added to the list of blocks to be processed L. 

[0394] The test and checking steps E51 3 to E51 9 are then reiterated over each of the blocks of the list L until the 

latter is empty 

[0395] In contrast, if. on conclusion of test E51 7, the detectability criterion is not satisfied on any of the child blocks 
of the parent block B, the tested block B is held in memory in a memory-storage step E509 as a reg.on available for 
ss inserting an information bit. 

[0396] In practice a step of updating the capacity E51 0 can be implemented as described previously, in such a way 
that, for each block of coefficients tested which is suitable for insertion of an information bit. the value of the capacity 
is increased by one unit. 
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[0397] Hence, on conclusion of these adaptive partitioning steps E511 to E519, an image insertion capacity can be 
determined, equal to Q bits, as can an appropriate partition of the image corresponding to this capacity O 
[0398] A watermark can then be formed, with length equal to Q bits, and can then be inserted, bit by bit. into the 
various regions of the partition. 

s [0399] Obviously, as before, the step E520 of insertion proper of the watermark can be coupled with the determination 
of the partition, in such a way that each block B available on conclusion of the detectability test E507, can actually be 
watermarked by the insertion of an information bit by modulation of the coefficients of this block B. 
[0400] The insertion step E520 is carried out in the same way as when partition is achieved by iterative combination. 
[0401] Moreover, in one enhanced embodiment of this technique of partition by iterative division, a division of a parent 

io block B into child blocks can be adopted on completion of the detectability test E507 only if it represents a gain in 
capacity, that is to say into a quadtree-mode division only if the detectability criterion is satisfied over at least two child 
blocks. Otherwise, it may be preferable to keep the parent block of larger size so as to guarantee better invisibility of 
the inserted watermark. 

[0402] Obviously, the partitioning techniques described previously are not in any way limiting. 
75 [0403] Notably, it is possible to use various iterative segmentation techniques by division/merging, to which there 
belong the quadtree block partitioning described previously, but also graph partitioning. In general terms, there are the 
so-called bottom-up methods and the top-down methods depending on whether the iterative partitioning means is 
merging or division. 

[0404] Hence it is possible, if appropriate, to combine two techniques of partitioning by division and combination. 
20 This is the case when an initial partition is chosen, (or example, intermediate between the two extremes which are, on 
the one hand, a single region and, on the other hand, a partition into regions of the smallest admissible size, allowing 
both the operations of division and of combination. 

[0405] It is also possible to implement the method of determining a partition in order to maximise the capacity of the 
digital image by using several different segmentation techniques and then choosing the partition which confers the 
25 greatest capacity Q for the digital image. 

[0406] It is also possible to determine the best partition possible, that is to say the one yielding a set of blocks of 
coefficients of largest possible size, for inserting a fixed-length message, less than the maximum capacity of the digital 
image to be watermarked. 

[0407] In all cases, the present invention makes it possible to generate and to insert a pseudo-random sequence 
50. once only over a set of coefficients which are available for carrying out the insertion, and, as a corollary, to apply only 
one single pre-distortion to the set of coefficients. 

[0408] It is thus possible rapidly to obtain a partition of the set of coefficients which makes it possible to insert a 
watermark with a fixed correct detection probability for an entire range of distortion lying between "no distortion' and 
"predetermined distortion" applied during the determination of the partition. Moreover, the simulation of maximum dis- 
ss tortion before insertion is a quasi-deterministic operation which confers an advantage with respect to a theoretical 
statistical calculation. 

[0409] A description will be given now of the insertion method according to the invention in a preferred embodiment. 
[0410] The watermarking signal is inserted after a transformation step E600, in a set of coefficients of the domain 
transformed into sub-bands. 

40 [0411] As illustrated in Figure 1 3, in a choosing step E601 , there is for example chosen the high-frequency sub-band 
HH,, with the first resolution level, corresponding to a high-pass filtering in the horizontal and vertical directions. 
[0412] There is thus a set of spatio-frequency coefficients, with a cardinal number equal to N, denoted for example 

X=< Xj, 1 <i<Nr. 

[0413] For an image I of size 512X512, the sub-band H^ is of size N = 256x256. 
45 [0414] Since the watermarking has to be imperceptible and indelible, and therefore difficult to locate and remove by 
pirating, it is chosen to insert a pseudo-random watermarking signal, spreading its spectrum in order to make this signal 
undetectable by spectral or statistical analysis. 

[0415] For example, a pseudo-random sequence w is considered, which follows a uniform law on the interval [-1,1] 
with w =\ w h 1 < i < Pr, the length P of the sequence being less than or equal to N. 
50 [0416] Naturally, any pseudo-random watermarking signal, with a known distribution and a nil mean, may be suitable. 
The most usual distributions for the watermarking signal w are, apart from the uniform distribution on [-1 : 1] mentioned 
above, the binary distribution^ -1,1 \ and the centred standardised Gaussian distribution N(0, 1). 
[0417] The modulation which it is wished to apply to the coefficients in order to insert the watermarking signal uses 
a linear model of the type: 
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X j ~ X j(i) + aW j ( Wlth 1 ~ j - P Where 
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[0418] X„ is a subset o, spectra, coe.tic.ents chosen from the set o, coefficients X. and o is a weight.ng coeff.c.ent. 
[0419] Here, the weighting coefficient^ is constant for any ■,. watermarking signal w is 

[0420] '"^-^ 

actually possible in the set of coefficients a wnusi r y h inc , udes a step E602 of calculating a maximum 

coefficients to be watermarked X. 

[0426] This calculation function used can be written: 



JS) 



where 
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is modulated, 

j} is strictly greater than 2, and 

E[lMrf]P is the mathematical expectation of the function lirvP. 



of the sub-band HH, is given by: 
so 6° 2 ct HH1 ( T) 



a v (N) : 



N 02 



^h^h^^ 

STtS« in order to meet a predetermined detectabi.ity criterion, the weighting coefficient a must be 
greater than a minimum value oq referred to as the detection amplitude. 
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[0435] At the decoder, the received image I*. which corresponds to the watermarked initial image, possibly made 
noisy because of certain distortions applied to the image, is first of all transformed by a transformation identical to that 
used during the insertion of the watermarking signal in order to find once again the set of coefficients which are mod- 
ulated. 

[0436] Thus a set X'HXj*. 1 <i<M is obtained, which represents the set of coefficients liable to have been modulated. 
By hypothesis, it is assumed that any distortion suffered by these coefficients can be modelled by an additive noise, 
decorrelated from the image signal itself n = \n h 1 < i < Nr.. 
[04371 Each coefficient can thus be written: 

X-* = X. + aw, + n.. with 1 < i < N. 



[0438] The detection of the inserted watermarking signal w consists of constructing a hypothesis test for replying to 
the question "has the signal w been inserted in X*?". 
is [0439] The advantage of constructing a hypothesis test is that its result makes it possible to give a theoretical con- 
fidence level to the decision since a hypothesis test is standardised with respect to known distributions. 
[0440] For a class of statistical tests based on correlation, it can be shown that this detection amplitude can be 
calculated from a calculation function of the type: 

20 n* — 2 

JbP - c 

where a£ is equal to the variance of the coefficients to be modulated, a 2 is equal to the variance of the additive 
25 noise and a, b and c are constants which depend on the distribution of the watermarking signal w and on the roquired 
detection probability ievei p. 

[0441] Thus this minimum value is determined according notably to the length P of the watermarking signal w, a 
detection probability level p, the variance of the coefficients to be modulated, the distribution of the watermarking signal 
and possibly the variance of an additive noise decorrelated from the coefficients modelling any distortions applied to 
30 the coefficients. 

[0442] Conversely, it is also possible to express the minimum detection length, denoted P°, as a function of the 
weighting coefficient a, for a fixed detection probability level: 

2 2 2 2 

35 P°(a) = — !-* f 

ba 

[0443] It would be possible to use this result to check the parameters which make it possible to insert a watermarking 
signal w which will be robust to any additive noise with an energy less than or equal to a n 2 . 

40 [0444] As illustrated on the curves in Figure 14, from the visual amplitude a v calculated as a function of the length 
P of the watermarking signal to be inserted, the insertion method includes, in general terms, a step of determining 
compatible values of the length P of the watermarking signal w, of the weighting coefficient a and of the probability 
level p so that the value of the weighting coefficient determined a is less than or equal to the maximum value or visual 
amplitude a v and greater than or equal to the minimum value or detection amplitude Oq for these determined values 

45 of the length P of the watermarking signal w and the detection probability level p. 

[0445] The weighting coefficient a, the length P of the watermarking signal w and the detection probability level p 
must thus be determined so that the following inequality is satisfied: 

so a 0 < a < a v 

[0446] As illustrated in Figure 14, in order to calculate the minimum value P° of the watermarking signal w for which 
the aforementioned inequality will be satisfied, the formula linking the length P of the watermarking signal w to the 
amplitude of the weighting coefficient a is used, with a = a^N) and with a detection probability level of p%. 
ss [0447] In this example a standardised statistical test is taken, based on the correlation between the watermarked 
signal X* and the inserted watermarking signal w, as proposed for example in the document entitled 'On Resolving 
Rightful Ownerships of Digital Images by Invisible Watermarks" by W. ZENG and B LIU, in Proc. ICIP 97, pages 
552-555, October 1997. 
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[0448] It can be shown that, for this test: 



[0449] A calculation step E604 thus makes it possible to calculate the three aforementioned constants a, b and c so 



10 that 



25 



30 



35 



45 
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a=t 0 JE(w 2 ), 

15 6=eW), 



and 

20 c=fo° 2 2- 



10450] The values E( **) and o 2 2 can be calculated theoretically as a function of the distribution of w. For the uniform 
distribution considered here, theoretical values arc E(v^)= 1/3 and <, 2 = _ 

ri^M? in this numerical example it is considered that any additive noise is not taken into account so that o - 0. 

Inlnesrat^ 
according to the following formula: 



40 with 



"'-Tit*' 

[0455] Then for a detection probability fixed at 99.95%, the function linking the length P o. the watermarking signal 
w to be inserted and the weighting coefficient a is then derived therefrom: 

so 14,44 X al+ 3,84a 2 

P= — 2 

0,11a 

,0456] A, « clcWion „ep E603, the nMmum len 9 lr, f is determined Iron MM abo.e toreieta to. a «eigM«9 
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the set of modulatable coefficients X, it is deduced therefrom that the insertion of a watermarking signal w of length P° 
is possible and meets the required imperceptibility and detection criteria 

[0459] In reality, it is possible to effect an insertion of a watermarking signal w which meets these criteria since the 
length P of this watermarking signal lies between the calculated value P° and the cardinal number N of the set of 
5 modulatable coefficients X. 

[0460] Then the insertion proper of a watermarking signal of length P is carried out by modulating a subset of coef- 
ficients with a cardinal number P according to the linear modulation level, with a weighting coefficient a equal to the 
calculated maximum value a v (N). 

[0461] The insertion method includes, in the usual fashion, a reverse spatio-frequency transformation step which 

10 makes it possible to obtain, at the output of the coder, the watermarked image I*. 

[0462] Where insertion is possible, it is also possible to calculate the number of theoretical bits which can be inserted. 
This is because, at the decoder, the decoding operation consists of detecting whether a pseudo-random sequence w 
has been inserted. The response being binary (yes/no), it can be considered that a pseudo-random sequence of length 
P makes it possible to insert and extract an information bit in the image. By repeating the insertion operation on several 

is subsets of spectral coefficients, it is thus possible to insert and extract several information bits. 

[0463] Consequently, since in this example there are a set of coefficients of cardinal number N and since only a 
number P° of coefficients are necessary to insert and detect an information bit, it is possible to derive therefrom the 
total number Q of information bits which can be inserted, referred to as the capacity of the image: Q = N/P°. 
[0464] Moreover, if at the comparison step E606 the calculated length P° is greater than the cardinal number N of 

■o the set of modulatable coefficients X, it is a priori deduced therefrom thai the insertion of a watermarking signal meeting 
the criteria of imperceptibility and detectability is not possible. 

[0465] In this case, the insertion method can if necessary include an additional step of calculating the detection 
probability level p' obtained for a minimum value ccq of the weighting coefficient equal to the maximum value <x v (N) 
calculated for a length P of the watermarking signal w equal to the cardinal number N of the set of modulatable coef- 



[0466] If this detection probability level, which wilt be less than that fixed initially at 99.95%, is judged satisfactory, 
the insertion of a watermarking signal of length equal to N can be effected, that is to say all the coefficients of the sub- 
band of frequency HH., are modulated. 

[0467] In a second embodiment of the invention, the insertion method can include a step of calculating a threshold 
value P* of the length of the watermarking signal w, as illustrated in Figure 14. 

[0468] This threshold value P* is determined so that the minimum value olq of the weighting coefficient is equal to 
the maximum value ct v of the weighting coefficient for a length P of the watermarking signal w equal to this threshold 
value P* and a predetermined detection probability level p. 

[0469] This calculation step makes it possible to determine, from formulae expressing the detection amplitude ct D 
and the visual amplitude a v notably as a function of the length P of the watermarking signal w, the minimum value P* 
of the insertion signal from which the weighting coefficient a can satisfy the inequality < a <a v . 
[0470] In Figure 14, this minimum value P* is given by the intersection of the two curves representing the variation 
in the visual amplitude a v and the variation in the detection amplitude as a function of the length P of the inserted 
watermarking signal w, for a predetermined detection probability level. 

[0471] Thus, in general terms, the insertion method according to the invention has the advantage of being able to 
control a priori and theoretically the detection probability and imperceptibility of the watermarking signal which it is 
sought to insert in an image. 

[0472] Naturally, many modifications could be made to the example embodiments described above without departing 
from the scope of the invention. 

[0473] Thus the insertion technique used could be applied to the raw digital image, without undergoing any spatio- 
frequency transformation prior to the modulation of the coefficients. 

[0474] In this case, the modulated coefficients are coefficients representing the digital image solely in the spatial 
domain. 

[0475] Moreover, the spatio-frequency transformation applied to the image can use analysis and synthesis filters 
other than those described above, or even be a transformation other than discrete Fourier transformation, discrete 
cosine transformation or Fourier-Mellin transformation. These transformations are used currently in conventional 
processings of digital images. 

[0476] In addition, the invention is not limited to the example embodiments described but also concerns any embod- 
iment where the values of the weighting coefficient a, of the length P of the watermarking signal and of the detection 
probability level p are determined in a correlated fashion so as to satisfy the inequality Oq < a <a v . 
[0477] The present invention also concerns a device for determining a partitioning and a decoding device which are 
adapted to implement the partitioning and decoding methods described previously. 

[0478] These devices for determining a partitioning and decoding can be used in a computer 50 as illustrated in 
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Figure 15 independently or in the same computer 50. 

[0479] The partitioning determination device comprises partitioning means 500, 502. 503 adapted to perform an 
adaptive partitioning according to a criterion of detectability of an information bit inserted on each region 
[0480] These partitioning means use on coding a maximum value a v of the weighting coefficient ensuring impercep- 
5 tibility of the watermarking signal for modulating all the coefficients representing said digital image, possibly with a 
modulation factor k, dependent on each coefficient, as described previously. 

[048i] This device for determining a partitioning also has means 500, 502, 503 of comparing the size P of each 
region of the partitioning with a minimum size P staf corresponding to minimum size of a statistically significant sample 
for the detection of a watermarking signal inserted in said region. 

10 [0482] The device also has means 500, 502, 503 for applying a pre<iistortion before calculating a delegability cri- 
terion. . , . 
[0483] The partitioning means comprise means 500, 502, 503 of modulating all the coefficients of each region by 
inserting the watermarking signal, means 500, 502, 503 of calculating a detectability amplitude from the detectabihty 
criterion and means 500, 502, 503 of validating the detectability criterion, adapted tocompare the detectability amplitude 

is with a predetermined threshold value Tc. 

[0484] The partitioning means 500, 502, 503 are adapted to effect a partitioning by iterative division of the numerical 
image and for each region at a rank in the said partrtioning, the validation means 500, 502, 503 are adapted to validate 
a partitioning of a rank higher in the region if and only it there exist at least two subregions of the region for which the 
detectability criterion is validated. 

20 [0485] This device also comprises means 500, 502, 503 of incrementing the value of a capacity Q of the digital image, 
adapted, for each region at a rank in the partitioning, to increment the capacity value Q when the partitioning of rank 
directly higher is not validated, the capacity Q being equal to the cardinal number of the distinct regions obtained by 
the partitioning on which the detectability criterion is validated. 

[0486] In another embodiment, this device for determining a partrtioning comprises incrementation means 500, 502, 
25 503 adapted for each region R at a rank in the partitioning, to increment the value of a capacity Q when the partitioning 
of rank directly higher is not validated, comparison means 500, 502, 503 adapted to compare the value of the capacity 
Q with a fixed capacity value Qs and reiteration means 500, 502, 503 adapted to reiterate the partitioning lor a region 
having the highest detectability amongst all the other regions to be processed when the capacity value Q is less than 
the fixed capacity value Qs. 

[0487] In a third embodiment ol the invention, the partitioning means 500, 502 : 503 are adapted to effect a partitioning 
by iterative merging of the digital image and there are adapted, for each region of the partitioning for which the de- 
tectability criterion is not validated, to merge the region with another region of the partitioning, and preferably with 
another region of the partitioning for which the detectability criterion is not validated. 

[0488] In an entirely similar fashion, the device for decoding a watermarking signal in a digital image comprises 
means 500 502, 503 of partitioning the digital image to be decoded into distinct regions, adapted to perform an adaptive 
partitioning as a function of a criterion of detectability of an information bit obtained by demodulating the coefficients 
on each region. 

[0489] It also has means 500, 502, 503 of comparing the size P of each region R of the partitioning with a minimum 
size P stat corresponding to the minimum size ol a statistically significant sample for the detection of a watermarking 
40 signal inserted in the region. 

[0490] It comprises means 500, 502, 503 of calculating a detection amplitude from the detectability criterion and 
means 500, 502, 503 of validating the detectability criterion, adapted to compare the detection amplitude with a pre- 
determined threshold value Td for decoding. 

[0491] This predetermined threshold value Td for decoding is less than the predetermined threshold value Tc used 
45 during the partitioning determination method described previously. 

[0492] In a preferred embodiment, the device for determining a partition includes means 500, 502, 503 tor checking 
a criterion of detectability, over each region, of an information bit inserted by means of -a pseudo-random sequence 
with average of zero by modulation of the coefficients of the region. 

[0493] It also includes means 500, 502, 503 for generating a centred pseudo-random sequence equal in size to the 
so cardinal number of the set of coefficients representative of an image, here equal to N. 

[0494] Means 500, 502, 503 of modulating the set of coefficients X by the centred pseudo-random sequence in order 
to insert the same information bit onto this set of coefficients are also incorporated into the device for determining a 
partition 11 as shown on figure 1. 

[0495] The latter may also includes means 500 : 502, 503 for applying a distortion D which makes it possible to 
ss provide for and guarantee the robustness of the inserted watermark in the case of certain types of distortion, typically 
a compression of the digital image in order to store it. 

[0496] These means of applying a distortion make it possible actually to simulate a predetermined distortion over 
the set of modulated coefficients. Here they comprise, whenever the coefficients of the set X are spectral coefficients. 
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means 500, 502, 503 for inverse spectral transformation of the set of coefficients modulated in order to obtain a wa- 
termarked image, means 500, 502, 503 for applying a distortion to this watermarked image, and means 500, 502. 503 
for spectral transformation of the watermarked image in order to recover the set of modulated coefficients after distor- 
tion. 

5 [0497] In this example, the inverse-spectral-transformation means are suitable for implementing multi-resolution 
spectral wavelet recomposition and the forward spectral transformation means are suitable for again reiterating the 
breakdown of the digital image I into wavelets. 

[0498] The device for determining a partition 11 , in a first way. also includes means for partitioning by iterative com- 
bination, on the basis of an initial segmentation of the set of coefficients X. 

10 [0499] For preference, this initial segmentation makes it possible to segment the set of coefficients X into a set of 
regions of a size which is statistically significant for validating the detectability criterion over these regions. 
[0500] The means of partitioning 500, 502, 503 by iterative combination are suitable, for each region of the partition 
for which the detectability criterion is not satisfied, for combining the region with another region of the partition. 
[0501] Obviously, the partitioning means 500, 502, 503 could also be suitable for forming this partition of the image 

is by iterative division of the digital image according to a second way of performing the preferred embodiment of the 
invention. 

[0502] This partition-determining device 11 makes it possible, as output, to obtain a partition into regions o1 the set 
of coefficients X, with it being possible to insert an information bit onto each of these regions with a reliable detection 
probability and robustness to a given distortion. The maximum number of regions which can be obtained by this adaptive 

20 partitioning of the set X corresponds to the practical insertion capacity of this set X. 

[0503] As illustrated in Figure 1, an insertion device 10, incorporating the device for determining a partition 11 de- 
scribed above, moreover comprises means 1 3 for reading a watermark to be inserted, containing several information 
bits. By way of an example which is not at all limiting, a binary watermark may be 1000111 , representing, for example, 
the coding of the name of the author of the image to be watermarked. 

25 [0504] The insertion means 1 2 proper may be conventional and make it possible, by modulation, to insert the various 
information bits of the watermark onto the various insertion carriers determined by the device for determining a partition 
11. 

[0505] According to the invention, the device for inserting a watermarking signal w comprises means 500, 502, 503 
of spatio-frequency transformation of an image I, and for example analysis filters associated with decimators by two 
30 adapted to effect a wavelet decomposition of an image I. It also has reverse spatio-frequency recomposition means 
500, 502, 503 for recomposing the image I after the insertion of the watermarking signal in the domain transformed 
into sub-bands. 

[0506] It also comprises means 500, 502, 503 of inserting a watermarking signal w of length P adapted to modulate 
a subset of coefficients of cardinal number P according to the linear model described previously using a weighting 
35 coefficient a. 

[0507] This insertion device also comprises means 500, 502, 503 of calculating the maximum value ct v of the weight- 
ing coefficient as a function of the length P of the watermarking signal w ensuring imperceptibility of the watermarking 
signal. These calculation means use the function a v (P, T) described above and- dependent on the length P of the 
watermarking signal w and the spectral transformation T used lor decomposing the image signal I. 

40 [0508] In general terms, this insertion device also has means 500, 502, 503 of determining compatible values of the 
length P of the watermarking signal w, the weighting coefficient a and the probability level p so that the value of the 
weighting coefficient determined is less than or equal tothe maximum value a v and greater than or equal to the minimum 
value Oq for the determined values of the length P of the watermarking signal w and the detection probability level p. 
[0509] The calculation means 500, 502, 503 are also adapted to determine a minimum value or detection amplitude 

45 oq using the calculation functioning described above of the type: 



[0510] In one embodiment of the invention, the insertion device also has means of calculating 500, 502, 503 the 
length P° of the watermarking signal wfor a minimum value oq of the weighting coefficient equal to a maximum value 
a v (N) calculated for a length P of the watermarking signals w equal to the cardinal number N of the set of modulatable 
coefficients ensuring imperceptibility of the watermarking signal w at the limit of perceptibility and for a predetermined 
detection probability level p; and means 500, 502, 503 of comparing this length P° calculated with the cardinal number 
N of the set of modulatable coefficients X. 

[0511] It also comprises means 500, 502, 503 of calculating a detection probability level p* obtained tor a minimum 
value Op of the weighting coefficient equal to the maximum value a v calculated for a length P of the watermarking 
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(0512] n f ^ 1 ^^^^^^t^ so that the minimum value <*o of the weighting coefficient is 
^eZ^gnt^^cent ot the length P of the watering signa, w eo.ua, to this 

5 r " ^^ESSSrS J^ST .^.^ SCO o, the computer SO, area. 
[ < Lmor! too befna adapted to store a program for determining a partitioning and/or .nsert.ng a watermarking 
SaTaT^ ^"t^SSSlinfl signa, and a random access memory 503 comprising registers adapted to 

,o rrr^^ 5 °' which can be connected to mefenx periphera,s d 

lot examp.e a ^SSS^S a mlrophone 511. by means o. an input/output card 506 in order to rece,ve and 
lOsTsnnTSgital camera 507 makes it possib.e notably to supply ,mages to be aulhenticaled by inserting a water- 
,s SieJ 9 ^computer 50 has a communication interlace 51 2 connected to a communication network 51 3 in order to 
5ST \Z:Z^^T^^ means, such as a hard d,sk 503. or is adapted tocooperate. by 

means of a disk drive 509. with removable document storage means such as diskettes 510. 

Sl8? These fixed or removable storage means can also contain the code of the drfferen, methods according to the 
invention which, once read by the microprocessor 500, will be stored on the hard disk 508 

invention whic^ ^ ^ ^ ^ ^ devjces tQ .^^^ ^ jnventjon cap be stored |n tne read only 

SSoT Afa second variant, the program can be received in order to be stored as described above by means of the 

^7 TheTo^SO also has a screen 504 for serving, for example, as an interface with an operator by moans 

« fc ^nT^ , in£ , 5^ execute the instructions relating to the ,mp.ementation of the invention On 
nSno uo ana methods relating to the invention stored with a non-vo.ati.e memory, for example the 

^r^^S'ZZn^ into the random access memory 503 (RAM) which wi.. then conta.n the exe- 
cutable code of the invention and the variables necessary for implementing the invent.on P11 „„ inn 
Sra ^random access memory 503 contains different registersforstoringthevariab.es ""^"S*^* 
oflne program a nd no tabfy a register .or storing the coefficients o. the regions R ^f^ZTl^^Z 
demodulated in order to verify the detectability criterion, the minimum value P Btat , and the threshold value Tc and Td 

[^Zr^^^ 503 also includes an * X, register for storing the coefficients of the regions 
Sare modulated or demodulated temporarily so as to check the detectability criterion and a reg.ster w for storing 
th P centred oseudo-random sequence generated in order to determine the partition of the set X. 
ZT fin^X T^ access memory 503 contains a register for storing the table of visual amplrtudes o U^. 
m a register for storing the spectral coefficients X, a register for calculating the visual amplitude « V (N) at the hm, of 
40 mperceS 

™r£Z^ l™ P yLe a,, as a function of the length P of the watermarking signa. w and a register «or storing the 

S^/I^^^SSIS. 501 affords communication between the different sub-elements o, the computer 50 or 
Sto a Threpresintation of the bus 501 is not limitative and notably the microprocessor 500 ,s able to .ndicate 
45 instructions to any sub-element directly or by means of another sub-element. 
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Method ol determining a partitioning into distinct regions of a digital image (I) in order to insert a watermarking 
slgnaUW by 3u.at?ng coefficients (X) of said regions (R) representing a digital image (I). charac,e«sed , , hat 
me paSoning into regions is effected by an adaptive partitioning as a function of a criterion of detectab.Uty of an 
information bit inserted on each region (R). 

Method of determining a partitioning according to Claim 1 , characterised in that the detectability criterion is calcu- 
latedCus^ 7* — value (« v ) of the weighting coefficient ensuring imperceptibly of the watermarking 
signal for modulating all the coefficients representing said digital image. 
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3. Method of determining a partitioning according to claim 2 characterised in that a maximum value (ct v ) can be 
calculated by: 

o. v (P, S, w) - — — - 

(E[lwlV 1 



where 

10 o^asefS) is a base value, dependent on the transformation used S and on the (base) sub-band considered for 

the insertion, of the maximum weighting coefficient ensuring imperceptibility when a single coefficient of said 
sub-band is modulated, 
P is strictly greater than 2, and 

E[M$] is the mathematical expectation of the function livll 3 . 

15 

4. Method of determining a partitioning according to one of Claims 2 or 3, characterised in that the detectability 
criterion is calculated by using, for each coefficient (Xj) to be modulated, a weighting coefficient in accordance with 
a law of the form a } = kj.a v , where k { is a modulation factor dependent on the coefficients situated close to the 
relevant coefficient (Xj) on said region and ay is equal to said maximum value of the weighting coefficient. 

20 

5. Method of determining a partitioning according to one of claims 1 to 4, characterised in that it comprises a step 
(E307) of applying a distortion before calculating the detectability criterion. 

6. Method of determining a partitioning according to one of Claims 1 to 5, characterised in that it includes a prior step 
25 (E105) of comparing the size (P) of each region (R) of the partitioning with a minimum size (P stat ) corresponding 

to the minimum size with a statistically significant sample for the detection of a watermarking signal inserted in 
said region. 

7. Method of determining a partitioning according to one of Claims 1 to 6, characterised in that, for each region (R) 
30 of the partitioning, all the coefficients of said region are modulated by inserting the watermarking signal (W), a 

detectability amplitude (T) is calculated from said detectability criterion and said detectability criterion is validated 
by comparing said detectability amplitude (T) with a predetermined threshold value (Tc). 

8. Method of determining a partitioning according to claim 7, characterised in that said detectability criterion is vali- 
35 dated it a minimum modulation length (P min (x)) calculated for ensuring detectability of the inserted watermarking 

signal (W) is less than or equal to the cardinal number of said region (R). 

9. Method of determining a partitioning according to claim 1 , characterised in that it comprises the following steps: 

40 - generation (E501 ) of a centred pseudo-random sequence (w) equal in size to the cardinal number (N) of a set 

of coefficients representative of at least a part of an image (I), formed from centred pseudo-random sub- 
sequences (w*); 

modulation (E502) of the said set of coefficients by the said centred pseudo-random sequence (w) in order to 
insert the same information bit on the said set of coefficients; and 
45 - checking (E506, E507, E514, E516), over each region, of a criterion of detectability of the information bit 

inserted by modulation. 

1 0. Method of determining a partitioning according to Claim 9. characterised in that it comprises a single step of applying 
a distortion (E503) to the said set of modulated coefficients, before the step of checking (E506, E507, E514, E516) 

so of a detectability criterion. 



11. Method ol determining a partitioning according to Claim 10, characterised in that the said set of coefficients is a 
set of spectral coefficients (X) originating from a spectral transformation of a digital image, and in that the said step 
(E503) of applying a distortion includes the following sub-steps: 

55 

inverse spectral transformation (E503a) of the said set of modulated coefficients in order to obtain a water- 
marked image; 

application (E503b) of a distortion to the said watermarked image; and 
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„,nn«ofli Claims 9 to 11 characterised in that the centred pseudo- 

said region. 

random samples generated, 
samples using a secret key. 

pseudo-random samples generated. 

mi n ir, ft n« of Claims 7 to 16 characterised in that a partitioning by 
17. Method ol determining a par1it.on.ng according to * 7 l ° ^ Qn at a rank of said partitioning, a 

said region for which said detectability criterion is validated 
«. Method of determining a partitioning according *»^$Z£^^ %Zt£S^££ 

srsrn^ 8 ^ 

by said partitioning on which the detectability criterion is validated. 

19 . Methodofdeterminingapari*^^ 

partitioning, the value of a capacity (Q) is incremented (E 11 2) when he J*™™ 9 partitioning 

when said capacity value is lower than said fixed capacity value. 

^ „ «^ ™ e of riaims 7 to 16 characterised in that a partitioning by 

20. Warned d = S SlaW* Ta^V .« toe^- ^on c saa P-Ui-*, ... «N* said 

ssr^^z - ^ - - — -*» - - 

a, Me^o^.-mininaapadaonin^^^ 

"7pSs» a "*on W wt*» said det.oabWy c™>"°« • ™» ™>"*><° & 

22 Ma.hod « .ns.dir* «,o a di g «a, «-g.. a — ""^J^^tT^ST S.tSS S 
said mg. in distinct ragioas forntina a paa««™>9 «■ ™* TliSi' no an adapts partition*., as a luno.,00 
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(P) of each region of the partitioning with a minimum size (P slal ) corresponding to the minimum size of a statistically 
significant sample for the detection of a watermarking signal inserted in said region. 

25 Decoding method according to one of Claims 23 or 24. characterised in that, for each region (R) of the partitioning, 
s ' a detection amplitude (T) is calculated (E203) from said detectability criterion and said detectability criterion is 

validated (E203) by comparing said detection amplitude (T) with a predetermined threshold value iTd) for decod.ng. 

26 Decoding method according to Claim 25, characterised in that the predetermined threshold value (Td) for decoding 
' is less than said predetermined threshold value (Tc) used during the partitioning determination method according 

10 to Claim 7. 

27 Device for determining a partitioning into distinct regions of a digital image for inserting a watermarking signal by 
' modulating coefficients of said regions representing a digital image, characterised in that it comprises partitioning 

means (500, 502, 503) adapted to effect an adaptive partitioning as a function of a criterion of detectability of an 
is information bit inserted on each region. 

28. Partitioning determination device according to claim 27, characterised in that it also includes means (500, 502, 
503) for applying a distortion before calculating the detectability criterion. 

20 29 Partitioning determination device according to one of claims 27 to 28. characterised in that it also has means (500, 
502 503) of comparing the size of each region of the partitioning with a minimum size corresponding to the minimum 
size of a statistically significant sample for the detection of a watermarking signal inserted in the said region. 

30 Partitioning determination device according to one of Claims 27 to 29. characterised in that the partitioning means 
25 (500 502 503) comprise means (500, 502, 503) of modulating all the coefficients of each region by .nsort.ng the 

watermarking signal, means (500, 502, 503) of calculating a detectability amplitude from said detectability criterion 
and means (500, 502, 503) of validating said detectability criterion, adapted to compare said detectab.l.ty amplitude 
with a predetermined threshold value. 

30 31. Partitioning determination device according to claim 27, characterised in that it includes: 

- means (500 502, 503) for generating a centred pseudo-random sequence (w) equal in size to the cardinal 
number (N) of a set of coefficients representative of at least a part of an image (I), formed by centred pseudo- 
random sub-sequences (w k ); 

35 - means (500. 502, 503) for modulation of the said set of coefficients by the said centred pseudo-random se- 

quence (w) in order to insert the same information bit on the said set of coefficients; and 

- means 500. 502, 503) for checking, over each region, a criterion of detectability of the information bit inserted 
by modulation. 

40 32. Partitioning determination device according to Claim 31 . characterised in that it comprises means (500, 502. 503) 
for applying a distortion to the said set of modulated coefficients. 

33 Partitioning determination device according to Claim 32. characterised in that the said set of coefficients is a set 
of spectral coefficients <X ( ) originating from a spectral transformation of a digital image, and in that the said means 

45 for applying a distortion comprise: 

- means (500, 502, 503) for inverse spectral transformation of the said set of modulated coefficients in order to 
obtain a watermarked image; 

- means (500 502, 503) for application of a distortion to the said watermarked image; and 

so - means (50o! 502, 503) for spectral transformation of the watermarked image in order to obtain a set of mod- 

ulated coefficients after distortion. 

34 Partitioning determination device according to one of Claims 30 to 33, characterised in that the partitioning means 
(500 502 503) are adapted to effect a partitioning by iterative division of the digital image and in that, for each 

ss region at a rank of said partitioning, the validation means are adapted to validate a partitioning of a directly higher 

rank of said region if and only if there exist at least two sub-regions of said regions for which said detectability 
criterion is validated. 
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47. Method of processing an image, characterised in that it comprises the following steps: 
inputting an image (E500): 

inserting a first watermarking information into said inputted image (E502): 
applying a distortion to the image obtained in said insertion step (E503a-c): 

detecting said first watermarking information from the image obtained in said distortion applying step 
(E504-E506); and 

determining an amount of a second watermarking information which can be inserted into said inputted image 
in accordance with the result of the detection obtained in said detection step (E510). 

48. Processing method according to Claim 47, characterised in that it futher comprises the step of displaying the 
amount of the watermarking information determined in said determining step. 

49. Image processing apparatus, characterised in that it comprises: 

75 

means (500, 502, 503) for inputting an image; 

means (500, 502, 503) for inserting a first watermarking information into said inputted image; 
means (500, 502, 503) for applying a distortion to the image obtained by said insertion means: 
means (500, 502, 503) for detecting said first watermarking information from the image obtained by said dis- 
20 lortion applying means: and 

means (500, 502, 503) for determining an amount of a second watermarking information which can be inserted 
into said inputted image in accordance with the result of the detection obtained by said detection means. 

50. Image processing apparatus according to Claim 49, characterised in that it further comprises the means (500, 502, 
2B 503) for displaying the amount of the watermarking information determined by said determining means. 

51. Method of inserting a watermarking signal (w) in a set of coefficients (X) representing a digital image (I), in which 
at least one subset of coefficients is modulated according to a linear model using a weighting coefficient (a), the 
weighting coefficient (a) being greater than a minimum value (oq) determined as a function notably of the length 

30 (P) of the watermarking signal (w) and a detection probability level (p), characterised in that it includes: 

a step (E602) of calculating a maximum value (a v (N)) of the weighting coefficients functional to the length (P) 
of the watermarking signal (w) ensuring imperceptibility of the watermarking signal. 

35 52. Watermarking signal insertion method according to Claim 51 , characterised in that it also includes: 

a step (E603) of determining compatible values of the length (P) of the watermarking signal (w) of the weighting 
coefficient (a) and of the detection probability level (p) so that the value of the weighting coefficient determined 
(a) is less than or equal to said maximum value (ay) and greater than or equal to said minimum value (ct D ) for 
40 said given values of the length (P) of the watermarking signal (w) and the detection probability level (p). 

53. Watermarking signal insertion method according to one of Claims 51 or 52, characterised in that, during the step 
(E602) of calculating the said maximum value (a v ), a function is used which is independent of said set of coefficients 
to be watermarked (X). 

45 

54. Watermarking signal insertion method according to one of Claims 51 to 53, characterised in that said coefficients 
are spatio-frequency coefficients obtained by spatio-frequency transformation of an image (I), the calculation of 
said maximum value (aj of the weighting coefficient depending on the type of transformation used (S). 

50 55. Watermarking signal insertion method according to Claim 54, characterised in that said calculation function used 
is of the type 

55 p 1/p (E[iwr]) 1/lJ 

where 
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56. 



57. 



a, (S) is a base value dependent on the transformation used (S) and on the (base) suM>and cons^ered 
fo^serlion of the max.mum weighting coefficients ensuring impercept.b.lity when a s.ngle coefhc.ent of s„d 

sub-band is modulated. 
P is strictly greater. than 2, and 

E[iwll ! ] is the mathematical expectation ol the function \Wf>. 

Watermarking signal insertion method according to one of Claims 51 to 55, characterised hMMnM»n 
value (an) is determined by using a calculation function also dependent on the variance (o *) of the coefficients 
^SmSZ distribution of the watermarking signal (w) and possibly the variance <«„*) «* an addit.ve no.se 
decorrelated from the coefficients modelling any distortions applied to said coefficients. 

Watermarking signal insertion method according to Claim 56, characterised in that said calculation function is of 
the type 

where o * is equal to the variance of the coefficients to be modulated, a* is equal to the variance of the 
additive noiae and a b and c are constants which depend on the distribution of the watermarking s.gnal (w) and 
the required detection probability level (p). 

58 Watermarking signal insertion method according to one of Claims 51 to 57, characterised in that the watermarking 
' signal (w) is a pseudo-random sequence with a predetermined distribution w.th a zero mean. 

59. Watermarking signal insertion method according to one of Claims 51 to 58, characterised in that it includes the 
following steps: 

. calculating (E602) a maximum value (o^N)) of the weighting coefficients for a length (P) of ^'^"^(j 
signal equal to the cardinal number (N) of the set of modulatable coefficients (X) ensuring imperceptibly of 
the watermarking signal (w) at the limit of perceptibility; Mainh ,, nn 

- calculating (E603) the length (P°) of the watermarking signal (w) for a minimal value (op) of the weighting 
Sdemequa. to said maximum value calculated (« V (N)) and for a predetermined detection probability level 

- Smearing (E606) said calculated length (P») with the cardinal number (N) of the set of modulatable coefficients 
(X). 
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Watermarking signal insertion method according to Claim 59, characterised in that, if at the comparison step i E606) 
%Z£££* length (P») is .ess than or equal to the cardinal number (N) of the set ° ^t^^TbTmS 
(X) there is inserted a watermarking signal (w) of length (P) at least equal to sad calculated length <P>)l by mod 
ulating a subset of coefficients with a cardinal number equal to said length (P) according to a linear mode. us.ng 
a weighting coefficient (a) equal to said calculated maximum value (a v (N)). 

61 Watermarking signal insertion method according to Claim 59, characterised in that, if at the °^|™ <^ 
me said length (>) is greater than the cardinal number (N) of the set of modulatable coefficients (XUhere is 
ScuEed the detection probability level (p-) obtained for a minimum va.ue (op) of the we.ght.ngj 

to said calculated maximum va.ue (a v (N)) for a length (P) of the watermarking signal (w) equal to the carnal 
number (N) of the set of modulatable coefficients (X). 

62 Watermarking signal insertion method according to one of Claims 5 1 to 58, characterised in that it 

of calculating a threshold va.ue (P*) of the .ength of the watermarking signal (w) determ.ned so that sa,d m.nimum 
value (oo) of the weighting coefficients is equal to said maximum value (« v ) of the weighting coefficient for a tengh 
(P)oi me Watermarking signal (w) equal to said threshold value (P') and a predetermined detection probab.hty 
level (p). 

63 Device for inserting a watermarking signal (w) in a set of coefficients representing a digital image (I) in which .at 
least one subset of coefficients is modulated according to a linear mode, us.ng a weighting coeff c,ent<« .the 
weighting coefficient being greater than a minimum value («o) determined as a function notably of the length <P) 
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of the watermarking signal (w) and a detection probability level (p), characterised in that it comprises means (500, 
502, 503) of calculating a maximum value (a v ) ol the weighting coefficient as a function of the length (P) of the 
watermarking signal ensuring imperceptibility of the watermarking signal (w). 

64. Watermarking signal insertion device according to Claim 63, characterised in that it also comprises means (500. 
' 502, 503) of determining compatible values of the length (P) of the watermarking signal, of the weighting coefficient 

(a) and of the detection probability level (p) so that the value of the weighting coefficient determined (a) is less 
than or equal to said maximum value (a v ) and greater than or equal to said minimum value (a D ) for said determined 
values of the length (P) of the watermarking signal and of the detection probability level (p). 

65. Watermarking signal insertion device according to one of Claims 63 or 64, characterised in that the means (500, 
502, 503) of calculating said maximum value (a v ) use a function independent of said set of coefficients to be 
watermarked. 

is 66. Watermarking signal insertion device according to one of Claims 63 to 65, characterised in that it comprises means 
(500, 502, 503) for the spatio-frequency transformation of an image so that the coefficients are spatio-frequency 
coefficients obtained by spatio-frequency transformation of an image, the calculation of said maximum value (a v ) 
of the weighting coefficient depending on the type of transformation used (S). 

20 67. Watermarking signal insertion device according to Claim 66, characterised in thai said calculation means (500, 
502, 503) use a function of the type 

a ba (S) 

a v (P, S, w) = — j ri7 

25 P 1/p (Ellivl P ]) 1/p 



where 

a base (S) is a base value dependent on the transformation used S and on the (base) sub-band considered for 
30 the insertion, of the maximum weighting coefficient ensuring perceptibility when a single coefficient of the said 

sub-band is modulated, 
P is strictly greater than 2, and 

E[lwlP] is the mathematical expectation of the function livlP. 

3S 68. Watermarking signal insertion device according to one of Claims 63 to 67, characterised in that said minimum 
value (a D ) is determined by using a calculation function also depending on the variance (o x 2 ) of the coefficients 
to be modulated, the distribution of the watermarking signal (w) and possibly the variance (c n 2 ) of an additive noise 
decorrelated from the coefficients modelling any distortions applied to said coefficients. 

40 69. Watermarking signal insertion device according to Claim 68, characterised in that said calculation function is of 
the type 



45 



a D (P) = 



a(Jo% 



JbP- c 



where c x 2 is equal to the variance of the coefficients to be modulated, o n 2 is equal to the variance of the 
additive noise and a, b and c are constants which depend on the distribution of the watermarking signal and the 
required detection probability level. 

so 

70. Watermarking signal insertion device according to one of Claims 63 to 69, characterised in that the watermarking 
signal (w) is a pseudo-random sequence with a predetermined distribution with a zero mean. 

71. Watermarking signal insertion device according to one of Claims 63 to 70, characterised in that it also comprises: 

55 

means (500, 502, 503) of calculating the length (P°) of the watermarking signal for a minimum value (c^) of 
the weighting coefficient equal to a maximum value (« V (N)) calculated for a length (P) of the watermarking 
signal equal to the cardinal number (N) of the set of modulatable coefficients ensuring imperceptibility of the 
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watermarking signal at the limit of perceptibility and for a predetermined detection probability level (p): and 
means (500, 502, 503) of comparing said calculated length (P°) with the cardinal number (N) of the set of 
modulatable coefficients (X). 

72. Watermarking signal insertion device according to Claim 71 , characterised in that it also comprises means (500. 
502, 503) of inserting a watermarking signal (w) with a length (P) at least equal to said calculated length (P°), 
adapted to modulate a subset of coefficients with a cardinal number equal to said length (P) according to a linear 
model using a weighting coefficient (a) equal to said calculated maximum value (a v (N), if said calculated length 
(P°) is less than or equal to the cardinal number (N) of the set of modulatable coefficients (X). 

73. Watermarking signal insertion device according to Claim 71, characterised in that it also comprises means (500, 
502, 503) of calculating a detection probability level (p') obtained for a minimum value (Oq) of the weighting coef- 
ficient equal to said maximum value (a v (N)) calculated for a length (P) of the watermarking signal equal to the 
cardinal number (N) of the set of modulatable coefficients. 



74. Watermarking signal insertion device according to one of Claims 63 to 70, characterised in that it comprises means 
(500, 502, 503) of calculating a threshold value (P*) of the length of the watermarking signal determined so that 
said minimum value (cxq) of the weighting coefficient is equal to said maximum value (a v ) of the weighting coefficient 
for a length (P) of the watermarking signal (w) equal to said threshold value (P*) and a predetermined detection 

20 probability level (p). 

75. Watermarking signal insertion device according to one of Claims 63 to 74, characterised in that it is incorporated 
in a microprocessor (500), a read only memory (502) storing a program for inserting a watermarking signal in an 
image and a random access memory (503) containing registers adapted to store variables modified during the 

25 running of said program. 

76. Computer, characterised in that it comprises means adapted to implement a partitioning determination method 
according to one o1 Claims 1 to 21 . 

30 77. Computer, characterised in that it comprises means adapted to implement the method of inserting a watermark 
according to Claim 22. 

78. Computer, characterised in that it comprises means adapted to implement the decoding method according to one 
of Claims 23 to 26. 

35 

79. Computer, characterised in that it comprises means adapted to implement a method of processing an image ac- 
cording to one of Claims 47 to 48. 

80. Computer, characterised in that it comprises means adapted to implement a method of inserting a watermarking 
40 signal according to one of Claims 51 to 62. 

81. Digital image processing apparatus, characterised in that it comprises means adapted to implement a partitioning 
determination method according to one of Claims 1 to 21 . 

is 82. Digital image processing apparatus, characterised in that it comprises means adapted to implement the method 
of inserting a watermark according to Claims 22. 

83. Digital image processing apparatus, characterised in that it comprises means adapted to implement the decoding 
method according to one of Claims 23 to 26. 

so 

84. Digital image processing apparatus, characterised in that it comprises means adapted to implement a method of 
processing an image according to one of Claims 47 or 48. 

85. Digital image processing apparatus, characterised in that it comprises means adapted to implement a method of 
55 inserting a watermarking signal according to one of Claims 51 to 62. 

86. Digital printer, characterised in that it comprises means adapted to implement a partitioning determination method 
according to one of Claims 1 to 21 
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87. Digital printer, characterised in that it comprises means adapted to implement the method of inserting a watermark 
according to Claim 22 

88. Digital printer characterised in that it comprises means adapted to implement the decoding method according to 
s one of Claims 23 to 26. 

89. Digital printer, characterised in that it comprises means adapted to implement a method of processing an image 
according to one of Claims 47 or 48. 

10 90. Digital printer, characterised in that it comprises means adapted to implement a method of inserting a watermarking 
signal according to one of Claims 51 to 62. 

91. Digital photographic apparatus, characterised in that it comprises means adapted to implement a partitioning de- 
termination method according to one of Claims 1 to 21 . 

15 

92. Digital photographic apparatus, characterised in that it comprises means adapted to implement the method of 
inserting a watermark according to Claim 22. 

93. Digital photographic apparatus, characterised in that it comprises means adapted to implement the decoding meth- 
20 od according to one of Claims 23 lo 26. 

94. Digital photographic apparatus, characterised in that it comprises means adapted to implement a method of 
processing an image according to one of Claims 47 or 48. 

25 95. Digital photographic apparatus, characterised in that it comprises means adapted to implement a method of in- 
serting a watermarking signal according to one of Claims 51 to 62. 

96. Digital camera, characterised in that it comprises means adapted to implement a partitioning determination method 
according to one o< Claims 1 to 21 . 

30 

97. Digital camera, characterised in that it comprises means adapted to implement the method of inserting a watermark 
according to Claim 22. 

98. Digital camera, characterised in that it comprises means adapted to implement the decoding method according to 
35 one of Claims 23 to 26. 

99. Digital camera, characterised in that it comprises means adapted to implement a method of processing an image 
according to one of Claims 47 or 48. 

40 100.Digital camera, characterised in that it comprises means adapted to implement a method of inserting a watermark- 
ing signal according to one of Claims 51 to 62. 

101 .Computer program product, characterised in that it is adapted to implement a partitioning determination method 
according to one ol Claims 1 to 21 . 

45 

102. Computer program product, characterised in that it is adapted to implement a method of inserting a watermark 
according to Claim 22. 

103. Computer program product, characterised in that it is adapted to implement a decoding method according to one 
so of Claims 23 to 26. 

104. Computer program product, characterised in that it is adapted to implement a method of processing an image 
according to one of Claims 47 or 48. 

ss 105.Computer program product, characterised in that it is adapted to implement a method of inserting a watermarking 
signal according to one of Claims 51 to 62. 

106.lnformation storage means readable by a microprocessor, characterised in that it stores a program implementing 
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a partitioning determination method according to one of Claims 1 to 21 

1 (^.Information storage means readable by a microprocessor, characterised in that it stores a program implementing 
a method of inserting a watermark according to Claim 22. 

5 108.lnformation storage means readable by a microprocessor, characterised in that it stores a program implementing 
a decoding method according to one of Claims 23 to 26. 

109.lnformation storage means readable by a microprocessor, characterised in that it stores a program implementing 
10 a method of processing an image according to one of Claims 47 or 48. 

110.lnformation storage means readable by a microprocessor, characterised in that it stores a program implementing 
a method of inserting a watermarking signal according to one of Claims 51 to 62. 
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